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A  numerical  procedure  is  developed  for  the  nonlinear 

analysis  of  the  structural  response  of  prestressed  concrete 

<i^.  members  subjected  to  fire.   The  fire  environment  is 

\::  simulated  by  the  time-tenperature  relationship  of  the 

standard  fire  test  (SFT)  .   The  procedure  con5>utes 
tenperature  distributions,  mid-span  deflection  and  moment 
capacity  of  a  member  supporting  its  design  service  loads 

'l   ,  ;•  '  •    throughout  the  simulated  fire  exposure  period. 

^' C'  .  The  formulation  of  the  numerical  model  and  the  solution 

of  the  equilibrium  equations  is  based  on  the  displacement 

■  *  'J 

finite  element  method.   The  member  cross  section  is 

idealized  into  triangular  elements  for  tenperature 

distribution  analysis  and  for  evaluation  of  variation  in 

material  properties  within  the  section .   Beam  elements  are 

used  to  model  the  member  for  both  deflection  and  moment 
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capacity  analysis.   The  transient  heat  transfer  mechanisms 
of  the  SFT  is  modeled  by  dividing  the  time  scale  into 
discrete  number  of  time  steps.   At  each  time  step 
equilibrium  equations  are  set  and  solved  iteratively  for 
tenperature  distribution  and  structural  response. 
1^  Incremental  deflections  within  a  time  step  are  added 

iv  cuimilatively  to  determine  the  total  deflection  at  the  end  of 

each  time  step. 

A  series  of  numerical  analyses  of  prestressed  concrete 
[  .         beams  and  slabs,  and  the  evaluation  of  the  influences  of 
tenperature-dependent  properties  of  concrete  and  steel  on 
teiT^serature  distribution,  deflection  and  moment  capacity  is 
presented.   Numerical  results  show  good  agreement  with 
corresponding  test  results .   Teirperature  distribution  is 
more  sensitive  to  thermal  capacity  properties  for  normal 
weight  concretes  and  to  thermal  conductivity  properties  for 
lightweight  concretes.   Deflection  is  mostly  influenced  by 
the  effects  of  tenperature  on  the  relaxation  properties  of 
prestressing  steel.   Thermal  properties  of  concrete 
influence  both  the  deflection  and  moment  capacity. 

Graphs  for  predicting  fire  endurance  of  prestressed 
concrete  beams  exposed  to  fire  have  been  developed  both  for 
normal  and  lightweight  concretes . 
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CHAPTER   1 
INTRODUCTION 


1.1.  General  Remarks 


The  determination  of  the  behavior  of  pres tressed 
concrete  members  subjected  to  high  tenperatures  has  been 
mostly  experimental .   The  experiments  are  based  on  the 
standard  fire  tests   (SFT) .   In  the  United  States  the  SFT 
procedures  are  defined  by  ASTM  E-119  [1] .   Results  from  the 
standard  fire  tests  have  been  used  to  develop  numerical 
procedures  for  fire  safety  designs. 

1.2.  The  Problem 


The  main  objective  of  the  Standard  Fire  Test  (SFT)  is 
to  predict  the  fire  endurance  of  building  con;>onents.   The 
pareuneters  measured  during  the  test  of  prestressed  concrete 
members  are  tenperature  distribution  within  the  member,  and 
the  mid-span  deflection  of  the  member  [2,3]  .   The  varicQsles 
considered  are  concrete  cover  and  type  of  concretes.   These 
variables  affect  the  tenperature  distribution  within  a 
member. 

Numerical  methods  developed  for  fire  safety  design  of 
prestressed  concrete  also  estimate  the  fire  endurance  time. 
The  parameter  used  is  the  moment  capacity  [3,4,5].   The 
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retained  moment  capacity  is  taken  as  a  function  of  the 
tenperature-dependent  tensile  strength  of  prestressing 
steel.   The  tenperature  of  steel  is  derived  from  tenperature 
distribution  curves  for  a  given  exposure  period  and  concrete 
cover.   However,  since  beams  are  exposed  to  fire  on  three 
sides,  these  curves  do  not  give  a  good  estimate  of 
temperature  of  steel  in  beams  [3] .   The  variables  considered 
are  the  same  as  those  of  fire  tests. 

Many  studies  show  that  almost  all  the  physical  and 
mechanical  properties  of  concrete  and  prestressing  steel 
undergo  considerable  changes  at  high  tenperatures 
[6,7,8,9,10].   The  procedures  of  the  standard  fire  test  does 
allow  for  the  evaluation  of  the  influence  of  these 
properties  in  the  measured  temperature  and  deflection 
behavior.   A  similar  problem  is  experienced  with  the 
existing  numerical  fire  design  methods . 

The  fire  tests  are  not  designed  to  measure  changes  in 
moment  capacity.   On  the  other  hand,  the  numerical  methods 
are  not  developed  for  predicting  changes  in  member 
deflection . 

Based  on  these  observations ,  it  is  paramount  to  develop 
a  method  that  will  address  all  the  above  short  comings .  The 
procedures  should  be  comprehensive  enough  to  incorporate  the 
changes  in  properties  of  concrete  and  steel,  and  predict  the 
history  of  temperature  distribution,  member  deflection,  and 
moment  capacity. 


1.3.  Research  Objectives 

The  objectives  of  this  study  are  as  follows: 

1 .  Develop  a  nxunerical  procedure  for  determining 
temperature  distribution  and  structural  response  of 
prestressed  concrete  members  exposed  to  fires.  The  procedure 
will  incorporate  the  nonlinearity  of  ten^aerature-dependent 
material  properties  of  concrete  and  prestressing  steel. 

2.  Use  the  procedure  developed  in  (1)  to  evaluate  the 
degree  of  influence  of  tenperature-dependent  properties  of 
concrete  and  prestressing  steel  on  the  confuted  teiq^erature 
distribution,  and  the  structural  response. 

The  histories  of  structure  deformation  and  strength 
capacity  are  considered  as  main  functions  for  defining 
structure  response.   The  numerical  model  developed  will 
evaluate  these  parameters  for  a  prestressed  concrete 
structure  exposed  to  a  simulated  fire  environment,  while 
supporting  its  designed  service  load. 

The  resulting  model  could  be  used  to 

(a)  estimate  teitperature  distribution  and  deformation 
status  of  a  loaded  prestressed  concrete  member  exposed  to 
fire.  This  application  will  reduce  dependence  on 
experiments , 

(b)  evaluate  strength  status  of  a  prestressed  concrete 
member  after  a  prolonged  exposure  to  fire  and 
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(c)  show  influence  of  tenperature-dependent  material 
properties  of  concrete  and  prestress  steel  on  temperature 
distribution  and  structural  behaviors . 

1.4.  Scope  of  the  Study 

The  tenperature-dependent  material  properties 
considered  in  this  study  are 

(a)  compressive  and  tensile  strength  of  concrete, 

(b)  stress-strain  relationship  of  both  concrete  and 
prestressing  steel, 

(c)  tensile  strength  of  prestressing  steel, 

(d)  thermal  strain  of  concrete  and  steel, 

(e)  shrinkage  strain  of  concrete  and 

(f )  accelerated  relaxation  of  prestressing  steel . 
Concrete  members  are  analyzed  on  the  assuitption  that 

the  moisture  content  is  within  the  equilibrium  relative 
humidity  of  70  percent  or  less.  This  state  siimilates  the 
common  conditions  attained  in  service  [11] . 

The  fire  environment  for  the  numerical  model  is 
simulated  by  the  time-teit^erature  relationship  defined  for 
the  SFT.   The  inposed  service  load  is  assumed  to  be  constant 
during  exposure  to  fire. 

Mathematical  models  of  tenperature-material  property 
relationships  are  derived  from  teitperature  test  results 
published  in  the  literature. 
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The  types  of  construction  materials  considered  are 
normal  weight  (silicon  aggregate  and  carbonate  aggregate) 
and  lightweight  (expanded  shale)  concretes  and  cold  drawn 
prestressing  steel  (wire  or  tendon) . 

1.5.  Problems  and  Solution  Schemes 


1.5.1.  Problem  Environment 

The  problem  involves  both  temperature  distribution  and 
structure  response  analyses.   The  main  conponents  of 
analysis  are 

(a)  tenperature  distribution  within  a  member, 

(b)  teitperature-dependent  structural  displacements  and 

(c)  tenperature-dependent  moment  capacity. 
1.5.1.1.  Teitperature  distribution  problem 

The  teitperature  distribution  problem  is  defined 
/^        graphically  in  Fig.  1.1.   Exposure  conditions  for  a  typical 
beam  element  are  shown  in  (a) .  The  problem  is  to  find  the 
teirperature  distribution  function  T(t)  under  the  transient 
heat  transfer  condition  of  the  SET  shown  in  (b) . 

Heat  is  transferred  from  the  source  to  the  element ' s 
exposed  surfaces  and  then  distributed  within  the  element. 
Heat  transfer  modes  involved  are 
-  (a)  radiation  and  convection  (Q  ^^  ^)  :    from  heat  source 

to  the  exposed  surfaces  of  an  element  and 


(b)  conduction  {Q^^J:   within  the  element. 
The  quantity  and  rate  of  heat  transferred  has  the 
following  dependencies : 


1  ^  h  c  ,  h  r 


(a)  Beam  section  showing  surfaces 
exposed  to  tire 


(b]    Temperature  models 


Fig.  1.1.  Tenperature  distribution  problem  environment 


(a)  Qhchr  is  a  function  of  heat  transfer  coefficients  of 
convection  (h ,,)  and  radiation  (h  ,.)  . 

(b)  Qk,pc  is  a  function  of  thermal  conductivity  (k)  and 
thermal  capacity  (  pc)  .  Higher  k  and  lower  qp  values  result 
in  higher  tenperatures  for  a  given  heat  transferred. 

The  above  variables  and  their  interrelation  have  to  be 
computed  continuously  in  the  process  of  solving  for  the 


temperature  distribution.   Detailed  procedures  for  the 
solution  of  tenperature  distribution  are  discussed  in 
Chapter  2. 
1.5.1.2.  Structure  response  problem 

Fig. 1.2  (a)  depicts  the  typical  system  configuration  to 
be  analyzed  for  structural  response.   The  system  is  under 
constant  service  load.   However,  due  to  effect  of 
tenperature  on  material  properties,  there  is  a  continuous 
redistribution  of  internal  stresses  within  an  element.   The 
core  of  the  problem  is  the  determination  of  changes  in  the 
internal  load  function  F,  and  the  response  function  r,  under 
the  constant  load  function  L,  as  shown  in  Fig. 1.2  (b) . 


LP 


L   C  <  on  3  t  ) 


P   C  T  .   e) 


L  =  applied  service   load 


P  =  function  of   prestresa  action 


stiffness   function   Cased 
on   properties   of   concrete 
and  prestressino  steel 


RA  =  reference  a» 


T  =   temperature  distribution   function 


e  =  stress-strain  function 


F     s      internal   load  function 


(t)   Loading  models  '      =      system  displncement   functi 

Fig.  1.2.  Structural  system  configuration 
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The  hypothetical  characteristics  of  the  main  parametric 
functions  of  the  analysis  are  illustrated  qualitatively  in 
Fig. 1.3.   Fig. 1.3  (a)  shows  the  decline  in  effective 
prestress  action  with  exposure  time  due  to  effects  shown  in 
Fig. 1.3  (b) ,  and  Fig. 1.3  (c)  shows  variations  of  stiffness 
parameter  S  and  displacement  r,  as  tenperature  T  increases. 
The  procedure  for  evaluating  these  parameters  is  discussed 
in  detail  in  Chapters  4  and  5 . 


L  .   P 


Ca)   Load   functions 


PC  n  3 


n  C  P ,   .  T.  ) 


service   load   function 


effective  prestress  action 


t  [  n  ) 


Cb)  prestress   load   functions 


f  CT  .   e  ,   0  .   g  ) 


=    function  of  prestress   loss 
due   to  relaxation  of   steel 

P  =   tension   in  prestress   steel 
=  f  (T  .      e  ) 

=    number   of   time  step 

=    function   of   prestress   loss 
due  to  geometric   changes 

T    =    Temperature   function 
t  =    exposure  time   function 


e   =   stress-strain   function 


f  (   S  .   R) 


f  C  T  ,   e  ) 


function  of   effectivt 
section  area 


Cc)  Stiffness,   displacement   functions 


Fig.  1.3.  Loading,  stiffness,  and  deformation  behavior 


9 

1.5.2.  Solution  Strategy 

The  entire  problem  involves  both  temperature 
distribution  and  structure  response  analysis .   The  problem 
is  formulated  and  solved  using  standard  finite  element 
techniques .   The  techniques  are  then  implemented  in  computer 
programs . 

The  prestressed  concrete  member  is  idealized  into 
finite  beam  elements  for  structural  response  analysis .   The 
cross  section  of  each  beam  element  is  modeled  into 
triangular  sub-elements .   The  cross  section  model  is 
utilized  for  tenperature  distribution  analysis,  evaluation 
of  the  variation  of  material  properties  across  the  section 
and  for  summation  of  the  overall  cross  section  properties. 
This  approach  is  illustrated  graphically  in  Fig.  1.4. 

The  time- temperature  relation  of  the  SFT  is  divided 
J  '  ,  into  a  number  of  thermal  time  segment.   The  time  boundaries 
of  the  segment  are  called  time  step  points.   Solutions  of 
equilibrium  equations  are  performed  at  these  time  steps  for 
both  tenperature  distribution  and  structure  response. 

The  solution  process  begins  with  solution  of  the 
temperature  distribution  problem.   Obtained  teir^serature 
values  are  used  to  confute  active  material  properties  for 
each  concrete  sub-element  and  prestressing  steel  segments, 
of  the  beam  element  section.   Based  on  the  variation  of 
tenperature  across  the  section  and  the  temperature-material 
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properties  relationship,  each  concrete  sub-element  and  steel 
segment  may  have  different  set  of  properties . 


y-'  ': 

\ 

\ 

frame     member 

[V\ 

-J^ 

\ 

\ 

\ 

/ 

\ 

\ 

*E^v- ' :  . 

/ 

'1 

^ 

member  modeled   into 
finite  beam  elements 


s  u  b  - 
e  I  e  m  e  n  t 


beam  element 
modeled  into 
subelements 


Fig.  1.4.  Finite  element  idealization  for  tenperature 
and  structure  response  analysis 


The  stiffness  parameters  of  concrete  sub-elements  and 
steel  segments  are  integrated  across  the  cross  section  of 
each  beam  element  to  determine  the  overall  element  stiffness 
characteristics .   The  strength  parameters  of  the 
prestressing  steel  are  added  across  the  section  to  determine 
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the  overall  effective  prestressing  effect  on  each  beam 
section.   The  confuted  functions  for  each  beam  element  are 
assembled  to  configure  the  structure  model.   The  solution 
process  is  then  continued  to  solve  for  structure  response. 

The  whole  procedure  is  inplemented  in  a  computer 
program  named  PRECET  (PRestressed  Concrete  at  Elevated 
Tenperature) .  The  program  consists  of  the  following 
sub-programs: 

(a)  TEDIAN  (TEnperature  Distribution  ANalysis) , 

(b)  STREAN  (STructure  REsponse  ANalysis)  and 
(b)  MAPRAN  (MAterial  PRoperty  ANalysis) . 
Discussion  on  the  development  of  the  conputer 

L  '^  .       inplementation  is  presented  in  Chapter  6.   The  program  is 

j  -  ■  ' 

\  ^-:  „       used  to  run  a  nuiiiber  of  exanples  which  have  been  tested  in 
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laboratories  according  to  Standard  Fire  Tests  (SFT) .   The 
numerical  procedure  developed  in  this  study  is  validated  by 
conparing  the  test  results  with  numerical  results .   The 
exanple  problems  investigated  and  the  conparison  of  results, 
are  presented  in  Chapter  7 . 

1.6.  Previous  Research 


Much  of  the  previous  research  related  to  structural 
response  of  prestressed  concrete  members  subjected  to 
building  fires  have  been  devoted  to  studying  the  effects  of 
high  temperature  on  the  properties  of  concretes  and  steels 
used  in  prestressed  concrete  structures  and  in  determining 
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the  tenperatures  within  a  concrete  member  during  the 
laboratory  standard  fire  tests  [3,5] , 

Behavior  of  concrete  and  prestressing  steel  under 
varying  teitperatures  have  been  studied  extensively  [5]  . 
Results  from  these  studies  show  common  qualitative 
behaviors,  though  may  differ  quantitatively  depending  on  the 
nature  of  materials  used  and  the  testing  procedures.   This 
section  presents  a  summary  of  findings  from  some  of  these 
studies.   More  information  regarding  material  properties 
behavior  will  be  given  in  Chapter  3. 

1.6.1.  Material  Property  Behavior . 

1.6.1.1.  Tensile  strength  of  prestressing  steel 

Day  et  al .  [3],  Abram  et  al.  [1]  and  Podolny  [10]  have 
studied  the  effect  of  temperature  on  the  tensile  strength  of 
prestressing  steel.   Results  show  that  the  strength  of 
prestressing  steel  decreases  by  more  than  50%  at  800  °F.   At 
1200°F  the  strength  is  down  to  5%  of  its  original  value. 

1.6.1.2.  Loss  of  prestress  force 

Abram  and  Cruz  [1]  observed  that  loss  of  prestressing 
amounted  to  15%  at  600  °F,   23%  at  700 °F  and  between  70%  and 
80%  at  1000 °F,  depending  on  the  initial  prestressing  force. 
Only  about  12%  of  the  loss  at  ten^jeratures  above  700  °F  is 
recovered  on  cooling.   The  test  was  performed  on  bare  wire, 
that  is,  not  embedded  in  concrete.   The  report  from  the 
study  did  not  isolate  the  parameters  that  caused  the  losses . 


/-•■. 


13 

1.6.1.3.  Cornpressive  strength  of  concrete 

Three  types  of  tests  have  been  used  to  determine  the 
relationship  between  the  con^sressive  strength  of  concrete 
and  teirperature  changes.   The  resulting  strength  categories 
are  as  follows : 

(a)  'unstressed'  strength  which  is  determined  by 
heating  the  specimen  with  no  superirtposed  loads.   When  the 
desired  test  teitperature  is  reached,  the  specimen  is  then 
tested  for  conpressive  strength. 

(b)  'stressed'  strength  which  is  determined  by  first 
loading  the  concrete  specimen  with  a  desired  test  load.   The 
stressed  specimen  is  heated  to  a  certain  tenperature  and 
then  tested  for  con^ressive  strength. 

(c)  'unstressed  residual'  strength  which  is  for 

]j>^   t '.  ■  '- 

[f         specimens  which  have  been  heated  to  any  desired  tenperature, 

».'  ■*"■  *  cooled  to  room  temperature  and  then  tested  for  conpressive 

strength. 

Results  from  compressive  strength  studies  by  Abrams  [9] 

and  Malhorta  [7]  show  that  most  concrete  types  retain  their 

near  full  strength  at  tenperatures  below  700   "F.  Siliceous 

aggregate  concretes  start  loosing  appreciable  amount  of 

strength  at  teitperatures  above  700   °F,  whereas  lightweight 

and  carbonate  concrete  losses  increase  at  tenperatures  above 

1000  °F.   However,  most  concretes  will  retain  only  about  20% 

of  their  room  temperature  strength  at  1600  °F.   It  was 

observed  that  type  of  aggregate  and  aggregate /cement  ratio 


J* 

r 
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have  the  most  significant  effect  on  the  strength  behavior. 
Specimen  tested  while  under  stress  showed  higher  values  of 
con^ressive  strength  than  for  unstressed  specimen. 
^  .         1.6.1.4.  Elastic  properties 
I  ^Sr,r  Elastic  properties  of  concrete  have  been  studied  by 

Cruz  [12]  and  Harada  [13],  among  others.   Results  from  Cruz 
f  show  that  the  elastic  modulus  of  concrete  drops  to  40-50%  at 

y.r  i       750°F,  and  below  30%  of  its  room  ten5>erature  value  at 

1200°F. 

Cruz  [12]  shows  that  the  modulus  of  elasticity  of 
^'  concrete  is  affected  primarily  by  the  same  factors  which 

^1  influence  the  con^jressive  strength.  The  aggregate  type  and 
'  concrete  strength  do  not  have  significant  effect.  However, 
as  for  conpressive  strength,  the  elasticities  from  stressed 
tests  show  higher  values  than  the  unstressed  ones. 
According  to  report  by  Harada  [13]  ,  original  room 
ten^serature  values  of  elastic  modulus  are  not  regained  upon 
cooling. 

The  effects  of  tenperature  on  the  elastic  modulus  of 
prestressing  steel  have  been  reported  in  Ref . [5] .   At  lower 
teir^jeratures ,  the  modulus  of  elasticity  decreases  at  a  slow 
rate  with  steel  teirperatures .   The  modulus  is  reduced  by 
about  6%  at  400 °F  and  by  approximately  20%  at  600  °F.   Above 
900°F,  the  modulus  of  elasticity  decreases  at  higher  rates 
to  almost  zero  at  1400  °F. 
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1.6.1.5.  Concrete  bond 


The  effect  of  terrperature  on  bond  strength  has  been 
studied  by  Morley  [14] ,  Morley  and  Royles  [15]  and 
Diederichs  and  Schneider  [16] .   Results  show  that  at  900  °F 
the  maximum  bond  strength  is  reduced  to  50%  and  to  less  than 
20%  at  1300 °F.   Effective  bond  between  concrete  and 
prestressing  steel  is  also  reduced  at  high  tenperature 
resulting  to  the  extension  of  development  length  [17,18] 
due  to  bond  slip.   Fagundo  and  Richardson  [19]  also  studied 
the  effects  of  tenperature  on  bond  slip  of  epoxy-coated 
ij         prestressing  strands.   Results  of  the  study  shows  a  bond 
j:-  „.  .,       slip  of  0.02  in.  when  the  strand  temperature  reached  of 
E         220°F.   The  strands  were  initially  stressed  to  75%  of  their 
ultimate  strength. 
'        1.6.1.6.  Deformation  and  thermal  properties 
Cl  Characteristics  of  nonmechanical  deformation 

properties,  such  as  concrete  shrinkage  and  thermal  strains 
of  concrete  and  prestressing  steel,  at  high  temperatures 
have  been  reported  in  various  studies.   Relationship  of 
thermal  properties  of  concrete  with  ten^jerature  have  also 
been  studied.   This  information  is  presented  in  Chapter  3. 
1.6.1.7.  Application  to  this  study 

The  models  of  tert^aerature-material  property 
relationships  used  in  this  study  have  been  derived  from  the 
results  published  in  the  literature.   Strength  properties. 
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nonmechanical  deformation  and  thermal  properties  of  concrete 
and  prestressing  steel  are  used  as  reported. 

Information  on  loss  of  prestress  is  not  used,  however, 
the  loss  of  prestress  is  derived  from  effects  of  concrete 
shrinkage,  thermal  strains  of  concrete  and  steel,  and 
relaxation  properties  of  prestressing  steel  at  high 
temperatures.   Modulus  of  elasticity  data  is  also  not 
applied  directly.   Elastic  properties  are  derived  from 
information  on  influence  of  temperature  on  stress-strain 
relation  of  concrete  and  prestressing  steel .   Due  to  lack  of 
sufficient  data  on  bond  slip  of  prestressing  steel  in 
relation  to  tenperature  changes,  bond  properties  are  not 
included  in  this  study.   Therefore  a  perfect  bond  between 
concrete  and  prestressing  steel  is  assumed  at  all 


*^\  '"--  tenperature  levels, 


1.6.2.  Structural  Response-Experimental 

Fire  tests  have  been  the  main  procedure  for  studying 
structural  response  of  prestressed  concrete  in  fire 
environment.   Most  of  the  test  data  reported  come  from 

(a)  Portland  Cement  Association  PCA) ,   Skokie,  and 

(b)  Underwriters'  Laboratories,  Chicago.   The  main  objective 
of  these  tests  was  to  determine  the  fire  endurance  of 
prestressed  concrete  members  while  supporting  their  full 
service  loads  [2,20,21,22]. 
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Fire  tests  of  building  con5>onents  are  conducted  under 
the  specifications  of  the  ASTM  Standard  Methods  of  Fire 
Tests  of  Building  Construction  and  Materials  [1] .  The  test 
san^le  is  placed  in  a  furnace  and  loaded  at  the  desired 
service  stage.   The  applied  load  is  the  maximum  permissible 
superinposed  load  as  required  or  permitted  by  recognized 
building  codes.   The  furnace  is  then  fired  by  natural  gas 
burners  from  each  side  of  the  sanple.   The  heating  of  the 
furnace  is  controlled  by  using  thermocouples  so  that  the 
tenperature  conforms  to  the  standard  time-tenperature  curve 
shown  in  Fig .1.5. 

Ten5>erature  distribution  within  the  specimen  is 
measured  through  thermocouples  located  at  various  points  in 
the  specimen.   Test  on  a  specimen  is  conducted  until  an 
arbitrary  selected  deflection  limit  or  a  defined  end  point 
is  reached.   Some  of  the  end  point  criteria  [1]  are 

(a)  Structural  end  point:  when  collapse  occurs  for  load 
bearing  specimens  while  supporting  the  applied  loading  and 

(b)  Heat  transmission  end  point:  when  the  tenperature 
of  the  unexposed  surface  of  floors,  roofs  or  walls  reach  an 
average  of  250  °F  or  a  maximum  of  325  °F  at  any  one  point. 

1.6.3.  Structural  Response-Numerical 

Research  areas  for  developing  numerical  procedures  have 
included  heat  transfer  and  structural  response  analysis 
[23]  .   Numerical  methods  developed  for  teirperature 
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distribution  include  the  use  of  finite  difference  techniques 
[13,24],  multi-degree  polynomials  [25],  and  finite  element 
[26,27] .   Finite  element  techniques  have  also  been  used  to 
develop  numerical  methods  for  structural  response  analysis 
of  steel  structure. 

Jeanes  [28]  developed  a  rational  method  and  conputer 
inplementation ,  FABUS  II,  for  determining  the  structural 
response  of  steel  buildings  subjected  to  fire.   Tenperature 
distribution  was  analyzed  using  FIRES-T3  [27]  .   Teir5>erature- 
dependent  properties  of  steel  considered  in  the  model  were 
the  yield  strength,  modulus  of  elasticity  and  coefficient  of 
expansion.   Structural  responses  studied  were  deflections, 
beeun  elongation  and  changes  in  stress  levels  across  beam 
sections . 

Anderberg  et  al.  [29]  also  developed  numerical 
procedures  for  predicting  behavior  of  steel  structures 
subjected  to  fire  and  iirplemented  in  a  FE  computer  program 
STEELFIRE.   Structure  responses  are  determined  by  the 
analysis  of  plane  steel  frames  subjected  to  in-plane 
loading.  Measured  temperature  values  were  part  of  the  input 
to  the  model.   Main  response  confuted  was  deflection. 

Other  developments  of  numerical  methods  for  response  of 
steel  structures  in  fire  include  those  of  Bresler  and  Iding 
[30]  and  Nakamura  et  al.  [31]. 

Numerical  methods  for  structural  response  of 
prestressed  concrete  structures  have  been  limited  to  the 
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determination  of  fire  endurance.   The  procedures  determine 
the  end  point  moment  capacity  of  a  prestressed  concrete 
member  subjected  to  elevated  ten^ieratures  [3,4,5,21]. 

The  residual  moment  capacity  is  confuted  using  the 
rectangular  beam  formula  given  in  the  concrete  building  code 
[32] .   The  basic  equation  is 


^ud  =  <\>Apsfps6  (d  -  ae/2)    =  (l)M„e 
a  = 


^8  ^pse 


0.85  fLh 


=  f. 


pi^ 


1  - 


ce^ 


0.5A_f, 


i.i: 


\ 


'ps    pu6 


bdfL 


c6 


where 

(|>  =  capacity  reduction  factor  for  flexure 
Ap,  =  the  cross  sectional  area  of  the  prestressing 

steel,  in. ^ 
fps  =  the  stress  in  the  prestressing  steel  at  the 

ultimate  load,  ksi 
fp„  =  the  ultimate  tensile  strength  of  the  prestressing 
steel,  ksi,  (teitperature  dependent) 
M„  =  nominal  moment  strength,  in.-k 
d  =  effective  depth,  in. 
a  =  the  depth  of  the  equivalent  rectangular  stress 

block  at  ultimate  load,  in. 
6  =  indicates  temperature  dependence 
Fire  safety  designs  aids  have  been  developed  based  on 
the  results  of  experimental  studies  on  ten^jerature 
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distribution  and  material  behavior.   The  aids  are  in  the 
form  of  graphs  for  teitperature  distribution  within  beams  and 
slabs,   and  for  temperature  and  strength  relationship  for 
both  concrete  and  prestressing  steel  [3,4,5,22,33,34]. 

Steel  temperature  within  the  concrete  at  certain 
distances  from  fire  exposed  surfaces  of  beams  or  slabs,  and 
at  various  times  of  exposure,  can  be  estimated  from  the 
teitperature  distribution  curves.   The  steel  strength,  f  p„, 
in  Eq.  (1.1),  is  then  confuted  according  to  the  relationship 
between  the  temperature  and  strength  of  prestressing  steel 
given  in  strength-tenperature  curves.   The  data  is  used  to 
compute  the  moment  capacity  at  various  exposure  times . 

Flexural  failure  is  assumed  to  occur  when  the  retained 
moment  capacity,  M^g,  is  reduced  to  the  level  of  the  applied 
moment  M.   This  condition  is  demonstrated  in  Fig.  1.5.   From 
Eq.  (1.1)  it  can  be  seen  that  the  fire  endurance  depends  on 
the  applied  loading,  and  on  the  strength- tenperature 
characteristics  of  the  prestressing  steel. 
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Design  service   load 


a)  Loaded  system 
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Service   load  moment,  M 


Moment   capacity  M 


b)  Moment  diagram  (before  exposure) 
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c)  Moment  diagram  (after   exposure) 


Fig.    1.5.   Moment  capacity  of  a  beam  exposed  to  fire 


.x« 


k-«  -'■■   -  i   ■- 

TF  ifcT  r-  ■  ,■■■ 


CHAPTER  2 
THERMAL   ANALYSIS 


2.1.  Introduction 


This  chapter  presents  the  development  of  the  procedure 
for  the  analysis  of  tenperature  distribution  history  of 
meiribers  in  a  fire  environment.   It  is  based  on  the  general 
approach  introduced  by  Wilson  [35,36]  and  Zienkiewicz  [37], 
and  simplified  for  the  needs  of  this  study. 

The  members  are  idealized  by  a  systems  of  two- 
dimensional  triangular  elements .   The  thermal  boundary 
conditions,  represented  by  the  time-teirperature  relationship 
of  the  Standard  Fire  Test  (SFT) ,  consists  of  convective  and 
radiative  heat  transfer  mechanisms .  Development  of  a 
conputer  program  (TEDIAN)  is  also  discussed. 

2.2.  Heat  Transfer  Models 


2.2.1.  Heat  Flow  Models 

The  general  transient  heat  transfer  equation  is 

a/,,ar\  ^   d i,dT\  _  „^ar  (2.1) 


-&m) '  i(*f )  -  "^f? 


where 

z,y  =  spatial  coordinates 
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T  =  tertperature  distribution  history 

t  =  time 

k  =  isotropic  conductivity  (temperature 

dependent) 
c  =  specific  heat  capacity  (temperature 

dependent) 
p  =  material  density  (teirperature  dependent) 
Initial  condition  is  specified  as 

T{x,y,C^o)    =  T^{x,y)  ^2. 2) 

where  Tq  is  the  specified  uniform  tenperature  for  every 
point  in  the  structure  before  analysis  begins.   In  this 
study  its  value  is  taken  as  68  °F.   Using  the  mathematical 
^1  5^  .  fc/   formulation  described  by  Segerlind  [38],  the  convective  and 
radiative  boundary  conditions  are  expressed  as 


;.>? 


where 

Ijj,  ly  =  directional  cosine  along  x-  and  y- 
coordinates ,  respectively 
hg  =  convective  heat  transfer  coefficient 
h^  =  radiative  heat  transfer  coefficient 
Tj,  =  surface  tenperature  of  concrete  member 
Tf  =  furnace  tert^ierature 
The  dependence  of  convective  heat  transfer  coefficient  on 
terrperature  has  minimal  effect  on  the  heat  transfer  process 
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*  "       of  the  standard  fire  test.   As  teinperature  increases  the 

t;         process  becomes  more  and  more  dependent  on  radiation.   For 
numerical  confutations  h  ^  is  assumed  constant  [39,40],  and 
approximated  as 

V  i2^  =  4.4   {Btu/hr.ft^°F)  ^^'^^ 

*  The  radiative  heat  transfer  coefficient  is  expressed  as 

h^  =  ae[{T^  +  460)2  +    (r^  +  460)2]  (r^  +  T^  +  920)         ^^'^^ 

where 

■ , «  o  =  Stef an-Boltzmann  constant 

-*  '■  '. 

=  0.119  *  10""  Btu/hr.in2.°R) 
H  •'  e  =  resultant  emissivity  of  the  flames,  combustion 

:^^        ,.,  '     .  gases  and  the  boundary  surface. 

:'"'■'  For  furnace  conditions  of  the  standard  fire  test,  the  value 

of  resultant  emissivity  is  taken  as  0.5  [39,40]. 

2.2.2.  Finite  Element  Formulations 


Based  on  variational  principles  [36,38,41],  Eqs.(2.1) 
and  (2.3)  are  expressed,  respectively,  as 


'  -  ilHf.! '  ifyf '  p^*^  -  ^p-*o* 


dV  (2.6) 


(2.7) 
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where 
■  '     [  -^  t(x,y,t)  =  unknown   tenperature  function 

'  *o  =  initial  ten^serature  field 

A    .  *s  =  tert^serature  field  of  surface  exposed  to 

heat  source  (unknown) 
*f  =  tenperature  field  of  medium  of  heat 
source  (known) 
The  above  functional  are  combined  so  that  the  boundary 
conditions  are  satisfied  automatically  after  minimization  of 
the  integral  I  [42] .   The  resulting  functional  is 


(2.8: 


2.2.3.  Matrix  Equations 

Thermal  matrices  are  formulated  using  standard 
displacement  finite  element  procedures.  Following  relations 
are  assumed: 


<|)(x,y,  t)  =  N{x,y)  T{t) 

where 

T(t)    =  nodal   ten5>erature  vector 
N(x,y)    =  interpolation   function 


(2.9) 


4 


TS 


<«^=(if) 


[D]    = 


k  0 
0  k 


[B] 


ox        dy 


(2.14) 


[H]    =    [B{x,y)]     (Tit)  ) 

Making  substitutions  of  Eqs.(2.10)  to  (2.13)  in 

^-   :   Eq.(2.8)  results  in 

^.»  .  •  ■    . 

:i      ;  (  I  =   {  \[T)-^[B]^[D]    [B]   [T)dV 

C^    .   ■'•     .  +  r  ^PC   {T)^[N]^[N]    {T)dV 

+    f  lh{T)''[N]''[N]   {T)  dS 
J  s  2, 

-   [  hTf[N]   {T)  dS  +   f  -hT^dS 
J  s  J  V  2 

where  h  represents  the  heat  transfer  coefficient, 
convective  or  radiative.   Eq. (2.14)  is  differentiated  with 
respect  to  tenperature  T,  and  equated  to  zero.   The 
minimized  function  is 
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(2.10) 

(2.11) 

(2.12) 
(2.13) 
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fi[B]^[D]    [B]dV  +  j  h[N]^[N]dsUT) 

+  f^pc[N]^[N]{{T)    -    {T,))dV  *2-^^^ 

-   r  hTf[N]  dS    =     0 

J  s 

For  an  integration  step     k,    tert^serature  rate  is  given  as 

where 

(TJ    =    (Tt)    =  initial   spatial   terrperature 
The  terms  in  Eq.    (2.15)    are  defined  as  follows: 

[K^]    =    (    [B]^[D]    [B]dV 

[JCJ    =    f  h^[N]^[N]  dS 
J  s 

[K,]    =    (  h,[N]^[N]dS 
J  s 

[C]    =   [  pc[N]^[N]  dV 
[i?J    =   f  h^Tf[N]  dS 

J  s 

[R^]    =   f  h^Tf[N]  dS 

J  s 


(2.17)      ;• 

(2.18) 

(2.19) 

(2.20) 

(2.21) 

i 
(2.22)         ■ 

■'I 


Eqs.(2.16)  to  (2.22)  are  lumped  into  a  matrix  form  of 
the  following  transient  heat  equation: 


^ 


^;;^ 


[c]  (r')    +   [K]  (T)    =   (J?) 

* .  N  where 

[K]    =    [KJ    +    [KJ    +    [K3]     (conductivity  matrix, 

teitperature-dependent) 
[R]    =    [Ril    +    [R2]     (thermal  load  matrix,    due  to 
pseudo  fire  boundary  conditions, 
temperature-dependent) 
[C]    =  thermal  capacitance  matrix,    teiqserature- 
f'  i   "  '  dependent . 
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(2.23) 


2.2.4.  Time  Integration  of  Tenperature 

Over  a  time  interval  (t+  k,t)  ,  the  rate  of  teitperature 
change,  (T'),  is  approximated  based  on  following  assunptions 
[42,43]: 

^t-aAt  =  ^(^t^At  -  Tt)  (2.24) 

rt.„At  =  (1  -  a)  r,  +  ar,,^t 

where  0  <  a  ^  1.   Then,  for  a  time  instant  (t+  atH  , 
Eq.(2.23)  is  rewritten  as 

f  ^J  t+aAt  (  ^  )  t+oAt    "•"     t-'^]  t+aAt^  ^^  t+aAt    ^     (■'^^t+aAt 

which  upon  substituting  in  Eq.  (2.24)  results  in 
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XE^^^t^-^t  -^  a[^]t.aAt 


(T) 


t+oAt 


:2.26; 


-  (1  -  a)  [i^]t.<xAt[r]t 


For  conditionally  stable  solution  of  Eq. (26)  in^licit 
algorithm  obtained  for  a  ''O  are  more  preferred.   Using  the 
algorithm  of  a =  1/2,  [35,42],  the  transient  heat  transfer 
equation  becomes 


-^[C],,^,  -H  [K],^ 


At 


At 


(2^)t.At=  (^)t.At  + 


;2.27) 


(i?)t+ 


^[c],.^,-[ic],.^,](r). 


Inplementing  and  trying  other  algorithms,  Eq.(2.27) 
required  more  confutation  effort  for  the  same  level  of 
accuracy  than  for  a=  1.  Thus,  the  form  of  equation  used  in 
this  development  is 


At 


[C-lt.At  +  [^]t.At  (r)t.At 


(2.28) 


-    (^)t.At  -^  -^^^C],^^,{T), 


The  more  convenient  form  for  solution  is 


[KC]   {r)t.At  =  (^c") 


(2.29) 


A 


where  the  effective  conductivity  matrix,  [KC] ,  is  given  by 
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[KC]    = 


^;^[C"],.^,  +  [^It.At 


(2.30; 


and  the  effective  thermal  load  vector,  [RC] ,  is 


[RC]  =   (i?),,^,  +  ^[c],.^,(r) 


;2.3i) 


n 


i 


re- 


2.3.  Evaluation  of  Thermal  Matrices  and  Load  Vector 


The  evaluation  of  thermal  matrices  of  Eq. (2.29)  is 
based  on  the  characteristics  defined  for  the  two  dimension 
triangular  element  presented  in  Fig.  2.1.  The  element  has 


■^    ■       one  tenperature  degree  of  freedom  per  node 


'>.,i 


Temperature   field 


>►  < 


spottAl   field 


Fig .2.1.  Plane  triangular  element 
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2.3.1.  Interpolation  Functions 

The  area  of  triangular  element  of  Fig.  2.2.  is  given  by 

The  three  shape  functions  are  defined  as  follows: 


■■,:-c. 


^i  =  ^[(^i:^ic-^icyi)  -  (yi-y^c)^  -  (^ic-^i)y]  '^•^^' 


^j  =  -2x[{''kyi-^iyk)  +  (yjc-yi)^  +  (^i-^jc)y] 


(2.34; 


iv, 


'^  ^  ^[(^i^i-^i^'i)  ^  (^i-^'i)^  ^  (^i-^i)y] 


;2.35; 


Substituting  Eqs.  (2.33) ,     (2.34)    and    (2.35)    in  Eq. 
(2.12)    we  get  matrix    [B] : 


^    ^        2A 


BI  BJ  BK 
CI  CJ  CK 


[2.361 


where , 


BI   =    (y^ 

-   Yx) 

BJ  =    (y. 

-   Yi) 

BK  =    (y. 

-  y.) 

CI  =  (x„  -  xj) 
CJ  =  (Xi  -  xj 
CK  =    (Xj    -   xj 


2.3.2.    Conductivity  Matrix 


i 


For  an  element  of  unit  thickness,  Eq.(2.17)  becomes 


F^i  =  i-A 


(BI^  +  CI^)  (BIBJ+CICJ)  (BIBK+CICK) 
(BJ^  +  CJ^)       (BJBK+  CJCK) 
SYM  (BK^  +  CK^) 


;2.37) 


Assuming  dS  =  1 . dL  and  using  area  coordinate  method  of 
integration  [44],  Eq.  (2.18)  takes  the  following  form: 
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F.]  =  ^ 


2 

1 

o" 

1 

2 

0 

.0 

0 

0. 

6 


0  0  0 
0  2  1 
0  12 


'2.38: 


where  i-j  denote  the  end  nodes  of  exposed  edge.   The  length 
L^j  is  defined  as 


Lij  =  slix^-x,)^  ^  iyj-Vi) 


;2.39) 


Eq. (2.19)  takes  the  same  form  as  Eq. (3.38) ,  replacing 
K   by  h^. 

2.3.3.  Thermal  Capacitance  Matrix 

Applying  the  area  coordinate  method  [42] ,  the 
capacitance  Eq.(2.20)  becomes 


[C]  =  -£^ 


12 


2   11 

2  1 

SYM        2 


(2.40; 


When  the  inplementation  of  Eq.(2.40)  was  used 
ten^erature  values  of  certain  nodes  diminished  below  zero 
before  rising.  The  reason  for  this  characteristic  is  not 
clear.   To  avoid  this  problem  an  alternative  approach  is 
utilized  in  this  study.   The  capacitance  matrix  is  con5>ut6d 
by  delineating  part  of  element  voliime  to  a  node.  For  element 
m  of  unit  thickness  the  volume  contribution  to  node  i  is 
given  by 


i 


m,  1  -^      m 


4A„p{r)  c{T) 


(2.41) 
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The  matrix  form  of  Eq. (2.41)    is 


I.     ••■        V      ;,->< 


[q  =  -£^ 


(2.42) 


The  approach  produced  results  with  the  expected 
characteristics . 

2.3.4.  Thermal  Load  Vector 


"   The  load  matrices  are  evaluated  using  line  integrals. 
For  an  element  of  unit  thickness  with  two  sides,  (i-j)  and 
(j-k) ,  exposed  to  fire,  the  evaluation  of  Eqs.  (2.21)  and 
(2.22)  result  to  a  con^jlete  thermal  load  vector  expressed  as 


iR)~-^^^I^ 


1 


^    KTfLj^ 


(2.431 


hrTfL, 


ij 


1 
\0. 


^     ^r^f^jk 


^0^ 

1 

VI, 


i 


2.4.  Solution  Algorithm 


At  time  (t+  At) ,  temperature  dependent  thermal  matrices 
[C] ,  and  [K] ,  of  Eq.(2.28),  are  not  known.  Thus,  for  each 
time  step  increment  It,    the  following  iterative  scheme  is 
used  to  obtain  an  acceptable  solution: 


'A 

: 
', 

1 


.*•>  ,;'"'"♦  u- 
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.-r\:.:-" 


..^ 


si^.-: 


At 


[q^^ 


(i-l) 
AC 


[^t. 


(i-l) 
At 


(^t+At  -   (■'?)t+At 


+  XtC^t'AV  (r)t 


where  the  superscript  (i-l)  denotes  iteration  number, 
initial  values  are  expressed  as 


;2.44) 


The 


^: 


[qt°At  =  [c]t 

r^(o)   _  r;^  (2.45; 

(i?)t°it  =(i?)t 

Convergence  of  solution  is  governed  by  both  the  system 
teiiperature  and  the  boundary  condition  criteria.   For  any 
two  successive  iterations  the  ten^jerature  convergence  is 
achieved  by 


|r7,(i)  _      rji{i-l] 

l|rp(i)  ,        rpd-l] 

■2|-'t+At  "^  -'t+At 


^  en 


(2.46) 


and  that  of  boundary  condition  by 


i^(i)  _h(i-i)| 


l|^,(i) 


^  e. 


(2.47; 
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CHAPTER  3  1 

MODELING  OF  TEMPERATURE  DEPENDENT  MATERIAL  PROPERTIES  ^ 

k 
3.1.  Introduction  ^jtj 


high  tenperature.   In  this  study  the  tenperature-dependent 
behavior  of  material  property  curves  is  estimated  by  a 
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Both  concrete  and  prestress  steel  experience  certain 
degree  of  degradation  when  exposed  to  high  temperatures 
[6,7,8,9].   It  is,  therefore,  inportant  to  understand  the 
behavior  of  these  properties  at  elevated  tenperatures . 

The  objective  of  this  chapter  is  to  present  the  ^.«;^ 

mathematical  relation  of  material  properties  to  .   ;\  :::S 

terrperature .   Most  of  the  models  are  derived  directly  from         j.| 
test  data  and  results  published  in  the  literature  and 
summarized  in  references  [45,5,46]. 

It  is  expected  that  the  data  used  and  the  models 
developed  represent  the  acceptable  behaviors  of  concrete  and 
|,'   ,.  i     prestress  steel  as  related  to  teitperature  and,  therefore, 
provide  reasonable  accurate  basis  for  evaluating  the 
structural  response  of  prestressed  concrete  elements   in 
fire . 

3.2.  Modeling  Procedures 

Most  properties  of  concrete  and  steel  are  affected  by         'f-i 


f' fl 


^,'' 
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series  of  discrete  points  connected  by  by  linear  segments  as 
shown  in  Fig.  3.1.   Then  the  parameter  value,  p,  at  a  given 
teirperature ,  T^,  is  confuted  by 


-■:S 


p(r^)  =p(r„)  +  (T^  -  T„)s„ 
5  ^  _p(r„.i)  -  PiTj 


T        -   T 

T     <.    T  ■    <    T 


(3.1) 


*■■' » 


••**■■ 


n 


f'-j 


P  (  T  )  A 


M 


Fig. 3.1.  Mathematical  Modeling  of  Material 
Properties  from  Test  Results 
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3.3.  Concrete 


'^■■'k    t 


3.3.1.  Thermal  Properties 

The  parameters  which  influence  the  characteristics  of 
heat  transfer  in  solids  are  the  thermal  conductivity  and 
thermal  capacity. 
3.3.1.1.  Thermal  conductivity 

Thermal  conductivity  of  concrete  measures  the  ability 
of  concrete  to  transfer  heat  by  conducting  heat  from 
particles  at  high  teitperature  to  those  at  low  teitperature . 
The  variation  of  thermal  conductivity  of  concrete  with 
tenperature  is  primarily  deteirmined  by  the  conductivity  of 
the  aggregates  used  [47,48].   High  tenperatures  increases 
the  disorder  of  particles  in  aggregates.   This  causes  the 
scattering  of  heat  in  all  directions,  hindering  effective 
transmission  of  heat  towards  low  teitperature  regions .   As  a 
consequence,  thermal  conductivity  of  concrete  decrease  with 
increase  in  tenperature . 

The  generic  model  relating  thermal  conductivity  of 
concrete  to  its  tenperature  is  expressed  as 

where 

Te  =  element  teirperature  (from  thermal  analysis) 
a,  b  =  constants  derived  from  experimental  data 
(according  to  Fig .  3.2.) 


"i 


The  variation  of  thermal  conductivities  with 
teirperature  for  the  three  types  of  concrete  is  shown 
graphical  Fig.  3.2. 
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200        400        600 


800       1000      1200      1400      1600      1800      2000 
Temperature   (deg.F) 


Fig.  3.2.  Conductivity  of  concrete  as  a 
function  of  tenperature 
(Source:  Ref.  [49]) 


3.3.1.2.  Thermal  capacity 

Thermal  capacity  of  a  material  refers  to  the  amount  of 
heat  required  to  raise  a  unit  volume  of  the  material  by  one 
degree.   It  is  a  product  of  material  density   (and  the 
specific  heat  of  the  material  c. 

The  character  of  the  cement  paste  and  water  in  the 
concrete  have  the  most  pronouncing  effect  on  the  behavior  of 
thermal  capacity  of  concrete  [47] .   In  the  case  of  water 
most  of  the  supplied  heat  energy  will  be  used  to  remove  pore 


1 

i 
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waters  in  temperatures  near  212  °F.   Crystalline  water  is 
removed  at  higher  tenperatures .   When  removal  of  water  takes 
place,  the  thermal  capacity  of  concrete  shows  substantially 
higher  values .   The  dependence  of  thermal  capacity  of 
concrete  on  terr^ierature  is  shown  graphically  in  Fig.  3.3. 
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Fig.  3.3.  Variation  of  thermal  capacity 
of  concrete  with  tenperature 
(Source:  Ref.  [49]) 


The  thermal  capacity  of  the  various  types  of  concrete 
considered  in  this  study  are  modeled  for  various  ranges  of 
tenperatures .   The  models  takes  the  following  form: 


oc[T)    =  a      +  b    T 


(3.3) 


1 
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3.3.2.  Mechanical  PropGrties 

Mechanical  properties  are  those  describing  the 
characteristic  of  material  under  stress.   Considered  in  this 
section  are  stress-strain  relation,  cort^sressive  strength  and 
modulus  of  rupture. 
3.3.2.1.  Stress-strain  relation 


High  tenperatures  increases  deformation  properties  of 
concrete.   Thus  there  is  a  continuous  variation  of  stress- 
strain  relation  with  ten^ierature  resulting  to  a  number  of 
strain-curves.  The  relationship  between  stress-strain  law 
and  terrperature  is  shown  in  Fig.  3.4  for  some  selected 
tenperature  ranges. 


-     -*.5 
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Concrele    Strain 

Fig.  3.4.  Stress-strain  relationship  as 
function  of  tenperature 
(Source:  Ref.  [49]) 
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For  each  ten5)eraturG  range,  modulus  of  elasticity  E  is 
confuted  for  a  given  range  m-n.   The  expression  is 

Ao^(r)  =  ojT)    -  ajT) 

!i  ^'■^;  ■  "/     '  ^^mn    =   ^n    -   ^m  (3.4) 


«:* 


l; 


■  5  ...' 


( ; 


~ntn 


-mi 


The  mechanical  stress  for  time  step  j+1,  temperature 
range  Tj^+i  and  mechanical  strain  €i  is  given  by 

^e^i  -  Ci    ^m  (3.5) 


The  stress-strain  curves  used  in  this  study  have  been  I 

reported  by  Lie  [49].   They  have  been  adjusted  to  account 
for  creep  by  moving  the  maxima  to  higher  strains  with  higher 
terr5>eratures .   In  addition,  the  curves  give  the  whole  i 

characteristic  of  stress-strain  relation  which  include  the 
post  peak  behavior. 

:  From  reported  data  [13],  aggregate  type  and  strength 

have  no  significant  effect  on  stress-strain  relation  at  high        ^ 
temperatures.   As  such  the  set  of  stress-strain  curves  J 

represents  both  normal  weight  and  lightweight  concretes .  | 

3.3.2.2.  Conpressive  strength 

'^}  Type  of  aggregates,  degree  of  loading  and  conditions  of        '^ 

loading-when  hot  or  cooled  are  eunong  the  inportant  factors  "^^ 

which  influence  the  behavior  of  concrete  strength  at 
elevated  temperatures. 


^  c: 


-■a 


'  < .  ■'  '■' ' 
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As  tenperature  of  concrete  increases  the  aggregates  and 

the  cement  matrix  expands.   At  the  same  time  pore  and 

crystalline  waters  are  evaporated,  which  causes  shrinkage  of 

concrete  matrix.   The  resulting  expansion  differentials 

cause  internal  cracking  of  concrete  which  results  to 

reduction  of  stiffness  of  concrete.   The  extent  of  this 

phenomenon  differs  considerably  with  the  type  of  aggregate  "^ 

used  in  concrete.  ] 

Concretes  with  silicious  aggregates  are  mostly  affected       j 

since  at  tenperatures  above  1000  °F  the  aggregates  also 

undergo  physical  changes  accornpanied  by  sudden  expansion  in 

volume  [50] .   Carbonate  aggregate  do  not  normally  undergo  ' 

physical  changes  during  heating,  hence  carbonate  aggregate 

concretes  are  free  from  severe  internal  cracks.   However,  at 

■'I 
very  high  temperatures  chemical  changes  take  place  when  lime        | 

i 

is  liberated  from  calcium  carbonate.   The  process  is  C| 

beneficial  in  retarding  temperature  rises,  but  during  5 

cooling  the  lime  combines  with  atmospheric  moisture  and 

expands  in  volume  causing  cracks  and  damages  to  the  concrete 

■1 
[50]  .   Lightweight  aggregates  normally  undergo  various  "i 

I 
heating  processes  during  manufacture,  hence  they  provide  a  '•• 

better  physical  con5>atibility  between  the  matrix  and  the 

aggregate  with  regard  to  deformation.     In  this  study  1 

experimental  data  used  to  derive  mathematical  models  are  for 

'stressed  tests'  in  which  concrete  is  tested  hot  while 

stressed  to  40  percent  of  its  room  ten^erature  strength. 


'  Tffli 
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These  results  are  meaningful  in  structural  response  analysis 
since  any  pres tressed  concrete  member  will  be  subjected  to 
some  stress  from  external  load,  prestress  or  both. 
According  to  Abrams  [9]  ,  the  coirpressive  strength  values  for 
specimens  stressed  at  25  to  55  percent  of  their  original 
compressive  strength  show  no  significant  variation  in 
behavior . 

The  graphic  models  of  the  variation  of  concrete 
strength  with  temperature  are  shown  in  Fig.  3.5, 
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Fig.  3.5.  Conpressive  strength  of  concrete 
versus  tenperature 
(Source:  Ref.  [5]) 


The  mathematical  expression  for  the  con^jressive  strength 
ratio  is  given  by  the  following  expression: 


I 


'^■■^r. 
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:3.6) 


CSR{T^)    =  C5J?(rJ  +  ^T^S^ 

where 

CSR  =  concrete  strength  ratio  (percentage  of  original) 
The  conpressive  strength  is  then  cornputed  by 

Fc(rj  =  C5i?(rj  Fc(rgg)      (3.7) 

3.3.2.3.  Tensile  strength 

Reported  test  results  of  tensile  strength  of  concrete 
at  high  temperature  are  very  sketch.   However,  these  reports 
'^'  indicate  that  the  deterioration  in  tension  is  greater  than 

',  "  that  in  conpression  [51]  .   For  lack  of  adequate  information, 
tensile  strength  of  concrete  is  taken  as  a  function  of 
compressive  strength.   The  relationship  is  expressed  as 

f^j.    =    7  .5  ^f^T  (  3  ^  g  ) 

fj..p   =5.6  ^fcT 

where  the  equations  are  for  normal  weight  and  lightweight 
concrete,  respectively  [32] . 

3.3.3.  Deformation  Properties 

Deformation  properties  considered  are  those  causing 
changes  in  the  dimension  of  a  member  due  to  non-mechanical 
means .   These  are  thermal  and  shrinkage  strains .   The 
thermal  strain  is  due  to  thermal  expansion  after  making 
allowance  for  drying  shrinkage.   This  is  only  an 
ass\imptionsince  according  to  Bazant  [52]  ,  it  is  difficult  to 


i 


.a 


..}  (.■ 


K  -i 


a^CT-)  =  — '^^^-^ 


^^th'^^mn)  (3.9) 


At 


Then  the  incremental  and  total  thermal  strain  at  temperature 
Ti  is  computed  from 

^th^'Ti)    =€,^(rj  +  Ae(7;,) 


■'ft 


45 
separate  thermal  strain  from  shrinkage  strain  under 
transient  tenperature  conditions.   Creep  strain,  which  is 
both  stress  and  time  dependent,  has  been  incorporated  in  the       -J 

derivation  of  stress-strain  relations.  -| 

J'. 

3.3.3.1.  Thermal  expansion  t| 

Thermal  expansion  of  concrete  increases  with  increasing        | 

teir5>erature .   The  variations  are  influenced  by  the  aggregate       '% 

;j 
type,  cement  content,  water  content  and  heating  rate  [51] .  J 

■A 
For  unstressed  specimen  the  increase  in  tenperature  causes  '• 

i 

the  increase  of  thermal  strain  to  a  point  of  disintegration. 
At  very  high  teirperatures  around  1100  °F  to  1500  °F  range  most 
concretes  indicate  no  expansion .  In  some  cases  the  concrete        '^ 
shrinks  due  to  chemical  or  physical  reactions  in  the 
aggregates . 

The  coefficient  of  expansion  is  derived  from  test 
results  of  thermal  strain  versus  tenperature.   The  relation 
is  shown  in  Fig.  3.6.   For  a  ten^serature  range  T  „,  T„,  the 
coefficient  of  thermal  expansion  is  given  by  the  following 
model: 
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Fig.  3.6.  Thermal  strain  of  concrete  as 
function  of  tenperature 
(Source:  Ref.  [5]) 


3.3.3.2.  Shrinkage  strain 

Shrinkage  strain  in  concrete  is  due  to  loss  of  moisture 
in  the  cement  paste  [52] .   The  amount  of  shrinkage  that 
occur  and  the  rate  of  its  occurrence  are  thus  dependence  on 
the  moisture  content  in  the  concrete  before  heating.     In 
this  study  the  concrete  members  are  assumed  to  be  at  normal 
service  condition  of  relative  htimidity  of  70  percent.   In 
addition  the  shrinkage  strain  considered  is  that  caused  by 
loss  of  free  water  from  the  pores  of  concrete.   Based  on 
this  criteria  all  the  shrinkage  is  assumed  to  take  place 
within  68  to  212  °F  teirperature  range. 
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The  incremental  shrinkage  strain  within  a  time  step  n 
is  confuted  by  the  following  models: 


Ae^^{r,)  =a{T,){ejT,)    -€,(rJ)At„ 


(3.11) 


a(r^)  =  0.001  + 


r^-6  8 

144 


(3.12) 


.if  7  J 


^a 


where 


e«(rj  =  0.0005  1  + 


(-4il^) 


(3.13) 


€s(Tn)  =  cumulative  shrinkage  strain  (at  beginning 
of  current  time  step) 
e.(T)  =  total  potential  shrinkage  due  to  loss  of 
free  water 
a(T)  =  shrinkage  rate  (in/in  per  hour) 
Aegh  =  incremental  shrinkage  strain 
At  =  duration  of  current  time  step  (hour) 
As  mentioned  earlier,  both  Eqs.  (3.12)  and  (3.13)  are  valid 
within  the  tenperature  range  of  68  °F  to  212 °F.   From  Eq. 
(3.13)  the  maximum  total  potential  shrinkage  at  212  °F  has 
the  value  of  0.001  in/in. 
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3.4.  PrestrGSsing  Steel 


Most  of  the  descriptions  of  property  categories  given 
in  for  concrete  apply  to  prestressing  steel  as  well. 


3.4.1.  Mechanical  Properties 

3.4.1.1.  Stress-strain  relationship 

For  a  given  range,  the  steel  modulus  can  be  obtained 
from  Eq.  (3.4)  after  substituting  for  the  appropriate 
variable.   The  steel  stress  is  then  given  by 


t*'   '  ■ 


Oi(r,e)  =  ajT^)   +  E^Ae„ 


(3.14) 


Fig.  3.7.  gives  the  graphical  representations  of  stress- 
strain  law  as  used  in  this  study. 
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Fig.  3.7.  Stress-strain  curves  for  PS 

at  various  temperature  levels 
(Source:  Ref.  [45]) 
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3.4.1.2.  Tensile  strength 
\     '  ,  The  tensile  strength  of  prestressing  steel  decreases 

'-•     with  the  increase  of  teitperature .   According  to  the  tests 
i^   ;,*  ?  '  done  by  Abrams  [6]  ,  the  size  of  the  strands  and  the  rate  of 

heating  have  no  significant  influence  on  the  tenperature- 

strength  behaviors.  Following  the  definitions  given  in  Fig. 

3.1,  the  tensile  strength  ratio  is  determined  by 


H 


'  •.  .<.. 


^TSR{T^)    =  TSRiTj    -  TSR{Tj 


At  ■  =  t-  -  t 


^mr,       ~ 


ATSRiT^) 


(3.15) 


TSRiT,)    =  TSRiTj    +  S^AT^ 


3.4.2.  Deformation  Properties 

3.4.2.1.  Thermal  strain 

Coefficient  of  expansion  of  prestress  steel  increases 
with  the  increase  of  ten5>erature .    It  is  expressed  in 
following  form: 

a(r)  =  (6.667  +  0.001235r)  (10"^)  ^  '  '^ ' 


s 


The  tensile  strength  is  obtained  from  V7I 

FPUiT^)    ^  TSR{T^)  FPUiT^)  (3.16) 

The  graphical  representation  of  the  teitperature- strength 
relationship  for  prestressing  steel  is  given  in  Fig.  3.8. 
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600  800  1000        1200        14-00         1600 

Temperature  (deg.F) 

Fig.  3.8.  Variation  of  prestress  steel 
strength  with  tenperature 
(Source:  Ref.  [5]) 


Variation  of  coefficient  of  thermal  expansion  of  prestress 
steel  with  ten^^erature  is  shown  in  Fig. 3. 9. 

Within  a  time  step  the  constant  value  of  coefficient  of 
thermal  expansion  used  in  confutations  is  given  by 


««  = 


The  corresponding  thermal  strain  is  expressed  as 


et/.(rj  =  e^  +  (^s^'^ni 


(3.18) 


(3.19; 


3.4.2.2.  Relaxation/creep  strain 

The  relaxation  property  of  prestressing  steel  can  be 
defined  as  the  decrease  in  stress  with  time  under  constant 
strain.   However,  most  of  the  reported  tests  have  measured 
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Fig.  3.9.  Expansion  coefficient  of  prestress 
steel  as  function  of  tenperature 
(Source:  Ref.  [45]) 


it  in  terms  of  rate  of  loss  of  stress  at  given  stress,  that 
is ,    creep . 

At  high  stresses  and  ten^eratures  creep  rate  become 
substantial.   According  to  tests  by  Yakovlev  et  al.  [53],  a 
prestressing  of  about  128,000  lb/in  ^  in  a  wire  was 
con^jletely  lost  when  heated  to  572  °F  in  40  minutes  at  a  rate 
equal  to  that  used  during  fire  resistance  tests  of  a 
prestressed  concrete  element. 

According  to  the  information  given  by  Harmathy  [54] , 
the  creep  characteristics  for  prestress  steel  are  described 
as  follows : 


II 


ir-s--;* 
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AH 
R 


=  54000 


Z  =  8.21  (10")  6°-°°°^" 


^cr  =  Ze 


RT 


(3.20) 


--> 


where 

AH/R  is  given  in  tert^jerature  units  (R) 

Z  =  Zener-Hollomon  parameter  (h  "^) 

T  =  absolute  tenperature  of  prestressing  steel  (R) 

o  =  stress  in  prestressing  steel  (lb/in  ^) 
ecr  =  creep  rate  (in/in. h) . 
Creep  rates  for  some  selected  stress  levels  are  presented  in 
Fig.  3.10. 
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Fig.  3.10.  Creep  rate  for  prestress  steel  as 
function  of  temperature  and  stress 
(Source:  Eq.  (3.20)) 
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3.5.  Remarks 

'  ? 

The  models  presented  in  this  chapter  have  been 

formulated  by  approximating  the  test  result  curves  by 

piecewise  linear  segments.   Similar  technique  was  used  by 

Franklin  [55]  to  model  stress-strain  relationship  of 

concrete  in  his  study  of  reinforced  concrete  frames .   This 

approach  was  used  in  this  study  because  of  lack  of 

mathematical  models  for  most  of  the  material  properties 

considered.   However,  generalized  enpirical  models  have  been 

used  whenever  possible. 


I 
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CHAPTER  4 
MODELING  OF  PRESTRESSED  CONCRETE  FRAME  MEMBER 


4.1.  Introduction 


'■i 


This  chapter  discusses  the  modeling  of  prestressed 
concrete  members.   The  formulation  is  based  on  direct 
stiffness  method  of  structural  analysis  using  two 
dimensional  beam  elements.   Using  the  idealization  of 
Fig.  1.4  the  stiffness  of  an  element  is  determined  by 
integrating  the  tenperature-dependent  subelement  properties . 

4.2.  Basic  Finite  Element  Equations  for  a  2D  Beam  Element 

Figure  4.1  presents  a  finite  planar  beam  element  with 
six  nodal  displacements,  r  ^  to  rg   These  displacement  J 

components  can  be  expressed  in  a  generic  form  as 

u  =  {u,v}''  (4.1) 

where  u  represents  displacement  cortponents  associated  with 

axial  action  in  the  x-axis  and  v  represents  both  translation        t, 

in  y-axis  and  rotation  due  to  flexural  action  about  the 

z-axis .   These  conponents  are  shown  in  Fig .  4.2. 

Relationship  between  the  generic  displacement   u  and  the 

nodal  displacements  [56]  is  given  by 
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u  ^  Nr 


(4.2) 


where  N  represents  the  displacement  shape  function,  and   r 
the  nodal  displacements  as  shown  in  Fig.  4.1. 


y 


t^ 


Fig. 4.1.    Displacement  cort^jonents  of 
a  2D  beam  element 
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Fig.  4.2.  Axial  and  flexural  elements 


Strain-displacement  relationship  is  expressed  as 


e  =  du 


4.3) 
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56 
where 

e  =  strain  vector 

d  =  linear  differential  operator 
By  Substituting  Eq.  (4.2)  in  Eq.  (4.3),  the  strain  equation 
becomes 

€  =  dNr  (4.4) 

resulting  to  the  definition  of  strain-displacement  matrix  B 

B-dN  (^-S) 


(4.6) 


The  corresponding  stress  is  given  by 

a  =  Ee 

=  EBr 

where 

E  =  stress-strain  matrix 
EB  =  stress-displacement  matrix 
The  relationship  between  the  nodal  displacements  and  the 
external  load  is  expressed  by 

where  K  is  the  element  stiffness  matrix  and   R  represents  the 
load  vector.   The  stiffness  matrix  is  formulated  as 


J   V 


B  ''EB  dv 


:4.8) 
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The  corresponding  load  vector  F,  due  to  strain  energy  of 
internal  stresses  in  the  element  is  confuted  from 


L 


F  =       B'^Eedv 


'A. 9] 


u  =  c^  +  c^x 


where  Ci  and  Cg  are  constants  corresponding  to  the  nodal 
displacements.   Expression  for  the  shape  function  is 


N  = 


-,  _  X  X 


;4.i3) 


L    L 
Strain-displacement  matrix  B  evaluates  to 


r 


'i4 


4.3.  Evaluation  of  Displacement  Functions  for  a  Beam  Element 

Based  on  the  actions  associated  with  the  nodal 

displacements  shown  in  Fig.  4.1,  the  beam  element  is  -  '   v  * 

deconposed  into:  axial  and  flexural  elements.   The  generic  J 

displacement  equation  for  the  axial  element  is  expressed  as  I 

where 

Corresponding  displacement  function  is  given  by  ^j 

(4.12) 


B=l[-ll]  (^-14)        , 

■■1 
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For  the  flexural  element,  displacement  relation  is  given  by 

The  corrponents  are  represented  by 

(4.16) 


r  = 


{r^,  z^,  r^,   r^}T  =  {v^,   0,,  v^,  e^}T 


where 


dx  .   (4.17) 

^        dx 

Based  on  plain  section  theorem,  translation  in  x  direction 
of  any  point  on  the  cross-section  of  the  element  is 
determined  from 

U=-yAX  (4.18) 

dx 

Substituting  Eqs.(4.18)  into  (4.3),  flexural  strain  becomes 

where  4)  represents  the  curvature.  Differential  operator 
takes  the  following  form: 


i 


d=  -y^L  (4.20)        -I 

dx^  ■'       ^ 
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The  displacement  function  is  expressed  as 


V  =    C^    +    C^X  +    CjX^    +  C4X- 


Shape  functions  are  then  given  by 


(4.21) 


i..- 


JT^  =  ^ 


1 


2x^   -  3x^L  +  L^ 

X^L  -  2x^L^   +  XL^ 

-2x^   +  3x^L 

x^L  -  x^L^ 


<■  (4.22) 


By  substituting  Eqs.(4.20)  and  (4.22)  in  Eq.(4.5),  the 
strain-displacement  matrix  B  becomes 


S^  =  - 


12x  -  6L 

/ 

SXL  -   4L2 

L^ 

-12j^  +  6L 

6xL  -   2L2 

(4.23) 


4.4.  Idealization  of  a  Prestressed  Concrete  Member 


A  prestressed  concrete  member  is  modeled  into  a  series 
of  beam  elements  and  tendon  segments  as  shown  in  the 
Fig.  4. 3  (a).   A  steel  segment  is  assiimed  to  be  straight,  has 
a  constant  tension  force  and  spans  the  beam  element.   The 
interaction  between  the  prestress  steel  and  concrete  is 
assumed  to  occur  at  the  ends  only. 

The  profile  and  the  location  of  a  steel  segment  in  a 
beam  element  is  defined  by  two  end  eccentricities  and  a 
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reference  plane  as  shown  in  Fig.  4.3(b).   It  is  determined 
by  the  expression  given  in  Eq. (4.24) . 


.-."i- '  * 


yp  =  T(^i  ^  ^j) 


(4.24) 


The  reference  plane  is  used  to  estimate  steel  segment 
ten^erature,  incremental  strains  and  contribution  to  the 
element  stiffness. 
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interaction  point  between 
and  steel  tendon 


Tendon  segmen  t  s 


a)  Finite  Element  Model  of  Beam  Specimen  and  Steel  Tendon 
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RA  =  reference  axis 
for  the  section 

RP  =  reference  plane 
for  steel 


b)  A  Beam  element  with  straight  tendon  segment 
and  a  constant  prestress  force 


Fig.  4.3.  FE  idealization  of  a  PC  member 
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4.5.  Evaluation  of  Element  Stiffness  Matrix 

The  strain  displaceinent  matrices  of  Eqs.  (4.14)  and 
(4.23)  are  modified  to  accommodate  all  six  displacement 
conponents  shown  in  Fig .  4.1.   The  matrices  become 


Sg  =  i.[-l  0  0  10  0] 


(4.25) 


for  the  axial  and 


-yBl  = 


y 


0 

12x  - 

6L 

6xL   - 

4L2 

0 

-12x 

1-  6L 

6xL  - 

2L2 

(4.26; 


for  the  flexural  element.  Eq.(4.8)  is  then  rewritten  as 


K  = 


=  /j(^a  -ySi,)''-E?(Sa  -  yB^)]dv 


or  as  the  following  summation: 


(4.27) 


^* 


4 


^-  ^aa    ^  ^ab^  ^ta    +  ^bb 


;4.28) 


where  the  conponents  are  expressed  as  follows : 


^aa=  [■ 

J    V 


B^EB^dv 


ai  =  /  -yB 


K 

Jc..  =  f  -yB 


aEBj^dv 


^ba 


i^EB  dv 


K^^  =   f  y^BlEB^dv 

J   V 


(4.29; 


^1 

.J 

■  .1 
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For  convenience  of  integration,    each  part  of  Eq.(4.29)    is 
decoitposed  in  two  integral   terms  as  follows: 

^ab  =   I  -yB^EB^dv  =  -  f  EydA  .  f  BlBj^dx  (4.30)  '■; 

K^i^  =  jy^BlEB^dv  =  J  Ey^dA.  J  BlB^dx  "    \  A'  ♦•      ; 

^ba        ^ab  ...  ^-        ,■•? 

The  variables  in  Eq.  (4.30)  have  following  properties: 
Matrix  B^  is  constant  as  shown  in  Eq.  (4.25).   Eq.  (4.26)  - 

shows  that  matrix  B^,  is  a  variable  of  x  only.   These  .  **;^ 

■   -'^ 
matrices  will  be  integrated  analytically.   The  stress-strain 

matrix  E  is  teitperature  dependent  and  a  function  of  both  x 

and  y.   In  this  study  each  subelement,  from  the  idealization 

of  the  beeun  element  in  thermal  analysis,  is  assumed  to  have         ^: 

a  uniform  temperature.   Thus  for  a  given  tenperature  E  is 

constant  for  each  subelement.   The  variable  A  is  the  section 

area  of  each  subelement  and  is  assumed  constant  along  the 

element  length  L.  Variable  y  gives  the  location  of  a 

subelement  or  the  projected  plane  of  the  prestressing  steel 

segment  with  respect  to  the  reference  axis  of  the  element. 

Based  on  the  assumptions  described  above  the  stiffness 

matrix  is  evaluated  term  by  term.   The  contribution  of  the 

prestressing  steel  to  the  element  stiffness  matrix  is 

obtained  by  considering  the  steel  segment  as  a  subelement 

located  at  its  projected  reference  plane  defined  in 


'■i 
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Fig,  4.3.  and  given  by  Eq.(4.24).  The  first  terms  of 
Eq.  (4.30)  are  integrated  niunerically  as  follows: 


psi 


EA   =  jj:dA  =  f;^  E,,A,,   .  g  E^^,A^ 
EQ  =  -JEydA  =  -^  E,,y,,A,,  -  "f^  E^,,y^^,A^^,  ^  ^  '  ^^  ^ 

^  ITl  i=l 

EI  =  JEy^dA  =   5^  E,,  (J,,  .y|,A,,)  -.  g  E^,,y^,,A^^, 


i-l 


1=1 


where 

^c/  ^ps  =  number  of  concrete  subelements  and  steel 
segments  respectively  within  an  element, 
Ap,  =  cross-section  area  of  prestressing  steel, 
loi  =  second  moment  of  area  of  a  subelement  about 
its  centroid  axis. 
Integration  of  the  second  terms  results  in 


L 


BiB^dx  = 


0 
0 

SYM 


0 
0 
0 


-1 
0 
0 

1 


(4.32) 


i  'J 
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L 


0    0 

1 

0    0 

0    0 

0 

0    0 

BlB^dx  = 

1 
L 

0    0 
0    0 

0 

-1 

0    0 
0    0 

0    0 

0 

0    0 

0    0 

0 

0    0 

-1 

0 
0 

1 

0 
0 


(4.33; 
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JBls.dx  .  A 


'00000         0 
12     6L     0    -12     6L 
4L2    0    -6L   2L2 

0   0    0 
SYM         12  -6L 
41,2 


(4.34) 


...  «.  / 


U" 


The  coitposite  element  stiffness  matrix  is  obtained  by   ' 
substituting  Eqs.  (4.32),  (4.33)  and  (4.34)  back  into       •  V'^  (f ; 
Eq.  (4.28)  .   The  resulting  matrix  is  given  by  Eq.  (4.35)  . 
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EQ 
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-12EI     -6  EI 
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-6  EI 
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L 

EQ 

L 

-6  EI 


L2 
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(4.35) 


4.6.  Incremental  Strains  in  Concrete  and  Prestressing  Steel 


4.6.1.  Concrete  Deformations 


Fig.  4.4  (a)  shows  the  incremental  displacements  of  the 
end  nodes  of  a  beam  element  after  two  successive  operations. 
These  displacements  are  with  respect  to  the  reference  axis 
(RA) .   The  centroid  axis  (CA)  changes  its  location 
continuously  due  to  changes  in  section  properties .   The 
incremental  displacements  are  transformed  to  relate  to  the 
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fb) 


i  ♦  I 


■1  c 


R 


(a  ) 


RA  =  reference  axi; 
CA  =  centroid  axis 


y   =  distance  of  CA 
form  R A 


1  =  iteration  number 


Fig.  4.4.  Incremental  nodal  displacements 


current  location  of  the  CA  before  corqsuting  the  strains. 
The  transformation  is  shown  in  Fig.  4.4  (b)  and  achieved  by 
following  expression: 


.CA 


(4.36) 


Using  the  transformed  values  the  effective  incremental 
displacements  translations  and  rotation  are  given  by 
Eqs.(4.37)  and  (4.38),  respectively. 


Ax  =  Ar^  -  Atj^ 
Ay  =  Atj  -  Ar2 


AG  =  Axg  -  Axj 
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(4.37) 

(4.38) 


Curvature  of  the  element  local  axis  is  expressed  as 

A(b  =  M  (4.39) 

Incremental  axial  strain  conponents  are  confuted  by 

^^a  =   4^  (4.40)       ^ 

for  axial  action  and 

for  bending  action.   L  „  is  the  original  element  length  and  y      -^  J 

is  the  distance  from  beam  element  centroid  axis  to  the 

centroid  of  the  subelement  in  considerations .   The  total  J 

incremental  strain  is  obtained  from  summation  of  Eqs.  (4.40) 

and  (4.41),  and  expressed  as 

Ae  =  A€,  +  A€^  (4-42)        :; 

Cumulative  incremental  strain  of  a  subelement  for  time  ^"^ 


step  i  and  iteration  j  is  given  by  Eq.  (4.43) 


J  _  ^i-i  .    ^    J  1 4.  431 


.(^ 
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4.6.2.  Prestressing  Steel  Strains 

Steel  strain  is  first  cort^iuted  for  the  projected  plane 
and  then  transfornved  to  the  correct  profile.  The  procedure 
uses  the  following  models : 

A€i  =  -ypA<|) 

Ae^-a  =  Ae^  +  Ae^  (4.44)  ■ 

Ae  =  Ae,_^  -    ^^^'^ 


cos  a 


The  cumulative  strain  for  a  steel  segment  is  similarly 
computed  from  Eq. (4.43) . 

4.7.  Load  Vector  due  to  Prestress  Force 

Element  nodal  forces  resulting  from  prestress  action 
are  shown  in  Fig .  4.5.   Following  the  definitions  given  in 
the  figure  the  length  of  the  steel  segment  and  its 
orientation  with  reference  to  the  x-axis  of  the  element  are 
computed  as  follows: 

Ax  =  XJ  -  XI 


Ay  =  e^. 

-  ©i 

'o  =  VAX^ 

+  Ay2 

cos  a  = 

AX 

sina  = 

AY 

■■•'1 


■r-< 

■   -J 


(4.45)        '' 

'i 
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where  XI  and  XJ  are  the  updated  x-coordinates  of  element 
nodes  i  and  j,  respectively,  and  L  „  is  the  original 
element  length. 


c:x^jL 


X  J 


R I  (  I  =  1  ,  .  .  ,  6 )  =  element  forces  due  to  prestress  action 


P  =  tension  in  the  tendon  segment 

P   ,Py   =  components  of  P  acting  on  concrete  at  interaction  point! 


i 


It' 


Fig. 4. 5.  Idealization  of  Prestress  Action 
into  Beam  Element  Nodal  Forces 


Conponents  of  the  prestress  force  along  local  axes  are 


P^  =   Pcosa 


Py  =  Psma 


:4.45) 


Translation  of  the  components  to  element  nodes  i  and  j 
result  in  the  following  end  moment  expressions : 
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R,   =  -Mj  =  P^ej 
The  element  load  vector  due  to  pres tress  action  is  given  by 

R     =  fp      p      -  P  e .  -P      -P      P  f^  \  ^   *  :  ■  '   (4.48) 


w 


4.8.  Remarks 


/■■:■'  f^ 


.■<»t^> 


The  formulation  of  the  beam  equations  in  this  chapter     *? ' 
is  based  on  the  standard  displacement  stiffness  method.   The     ''' \% 
evaluation  of  the  stiffness  matrix  is  performed  term  by  term 
across  the  section  element.   This  approach  is  applied  in 
order  to  incorporate  the  variations  of  stiffness  properties 
across  the  section  due  to  the  state  of  temperature 
distribution  and  stress-strain  relationship  within  different 
parts  of  the  section. 


li. 


CHAPTER  5 
SOLUTION  METHOD  FOR  STRUCTURAL  RESPONSE 


5.1.  Introduction 


This  chapter  presents  the  technique  for  solving 
structural  responses  of  prestressed  concrete  members 
subjected  to  high  tenperatures .   The  basic  equilibrium 
equations  to  be  solved  are  given  as  follows : 

(a)  At  the  beginning  of  a  solution  cycle 

(a)    Kr  =  R 
ib)  JTAr  =  A12 


:5.i) 


where  the  stiffness  function  K  and  the  applied  load  function 
R  or  AR  are  known.   The  task  is  to  solve  for  the 
displacement  function   r  or  Ar. 

(b)  At  the  end  of  a  solution  cycle 

(a)    Kr  =  F 

(5.2) 
{h)  KAr  ^  AF 

i 

.'  a'' 

In  this  case  the  stiffness  function   K  is  an  update  based  on 
the  displacement  function  r  or  Ar  computed  in  Eq.  (5.1) . 
The  internal  load  function  F  is  the  unknown  to  be 
determined. 

Techniques  for  evaluating  the  stiffness  function   K  have 
been  discussed  in  Chapter  4.  The  evaluation  of  the  load 
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functions  R  and  F  requires  the  continuous  determination  of 
the  following  parameters  in  all  solution  cycles: 

(a)  Induction  forces  and 

(b)  effective  prestress  action 
These  parameters  are  discussed  below. 

5.2.  Induction  Forces 


Induction  forces  are  defined  here  as  those  forces  which 
are  caused  by  stress  induced  in  the  system  as  a  result  of 
the  effects  of  tenperature  on  properties  of  concrete  and 
prestressing  steel.   The  induction  force  vector   AG  is 
given  by 


AG  = 


-A 
0 

AaJyA 
A 
0 
y-yA] 

Ao  =  Aa^  +  Aa^  +  Aa^ 


where  Ao  is  the  vector  of  the  induced  stress.   The  main 
contributing  factors  to  the  induced  stress  function  are 

(a)  thermal  and  shrinkage  strain  of  concrete  (  Ao')  , 

(b)  strength  degradation  (  Ao'^)  and 

(c)  shift  in  stress-strain  relation  (  Ao*)  . 


(5.3) 


M 
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5.2.1.  Thermal  and  Shrinkage  Strain 

Stresses  induced  to  the  system  due  thermal  and 
shrinkage  strains  are  the  equivalent  mechanical  stresses 
required  to  effect  equal  strains  if  the  system  was  free  to 
undergo  the  deformation.   Thermal  strains  of  concrete  are 
cort^juted  using  Eq.  (3.10)  and  shrinkage  strains  by 
Eq.  (3.11).   For  each  section  subelement  the  induced  thermal 
and  shrinkage  stresses  at  time  step  i  are  given  by 

where 

AOj^  =  induced  stresses 
ACi  =  thermal  and  shrinkage  strain 
Ei  =  valid  modulus 
Since  the  teit?>erature  of  a  subelement  does  not  change 
within  a  time  step  stresses  due  to  thermal  and  shrinkage 
strains  remain  constant  in  all  operations  within  a  time 
step.   They  are  computed  at  the  beginning  of  the 
time  step  only. 

5.2.2.  Strength  Degradation 


I'j. 


^ 


Based  on  the  material  property- teir^jerature  relationship 
both  concrete  and  prestressing  steel  degrade  as  temperature 

increases.   Degradation  is  defined  here  as  the  deterioration        -;, 

■/  ■  ;';''!| 

in  strength  and  stiffness  properties.   A  concrete  subelement        '^ 
is  considered  degraded  when 
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(a)  crushed:   o^  ^  f  c(Ti)  determined  by  Eq.  (3.7) 

(b)  cracked:   Oi  >ft(Ti)  determined  by  Eq.  (3.8) 

(c)  Eo  :£  0  as  determined  by  Eq.  (3.4). 
Criteria  for  prestress  steel  degradation  are 

(a)  tension  failure:  f  p.  ^  fpu(Ti)  determined  by 
Eq.(3.16) 

(b)  E,  ^  0  as  determined  by  Eq.  (3.4). 

Concrete  subelements  and  steel  segments  flagged  as 
degraded  are  excluded  in  the  corr^jutation  of  section 
properties.   The  stress  which  was  previously  carried  by  a 
degraded  concrete  subelement  is  redistributed  into  the 
system.   The  above  parameters  controlling  the  material 
degradation  are  both  ten5>erature  and  strain  dependent. 
Since  total  strain  may  change  within  a  time  step,  these 
induced  stresses  are  corrputed  for  each  iteration.   The 
stresses  induced  into  the  system  due  to  degradation  are 
expressed  as 

Aaf.t^  =   A  oil  (5.5; 

where  i  is  the  time  step  and  j  is  the  iteration  number. 
Superscript  D  is  an  indicator  for  degradation  induced 
function.   For  prestressing  steel,  the  degraded  steel 
segments  do  not  contribute  to  the  prestress  action  of  the 
containing  beam  element. 
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5.2.3.  Shift  in  stress 


At  the  beginning  of  a  time  step  the  active  stress- 
strain  curve  may  change  due  to  new  teitperature  range.   The  * 
resulting  effect  is  that  a  concrete  subelement  will  support 
less  stress  though  the  strain  remains  the  same.   This 
phenomenon  is  illustrated  in  Fig.  5.1.   The  stress  induced 
in  the  system  due  to  this  concrete  behavior  is  given  by 


AaU  =   ai  -  aU        ''    ^  ■  "  '        -'^    (5.6; 
where  S  indicates  source  as  stress  shift,  and  j  is  the  last 
iteration  of  the  previous  time  step. 


6  6(1,) 


A6i., 


Fig.  5.1.  Shift  of  stress-strain  curve 
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5.3.  Effective  Prestress  Action 

Prestress  action  in  a  beam  element  is  dependent  on  the 
effective  tension  in  the  contained  steel  segments.   The 
tension  in  a  steel  segment  is  a  function  of  total  steel 
strain  for  a  particular  teitperature  range.   Since  strain 
changes  with  each  iteration,  steel  tension  has  to  be  updated 
continuously  within  each  time  step. 

The  theoretical  steel  tension  is  subject  to  reduction 
due  to  partial  loss  of  stress.   The  function  considered  in 
this  study  as  the  main  sources  of  loss  of  prestress  are 
stress  relaxation,  thermal  and  shrinkage  strains  of  concrete 
and  thermal  strain  of  steel.   The  strains  caused  by  these 
parameters  are  assumed  to  remain  constant  within  a  time  step 
hence  are  determined  at  beginning  of  each  time  step.   Though 
steel  modulus  and  stress  may  change  within  a  time  step,  for 
sin5>licity  the  partial  loss  of  stress  is  held  constant 
within  a  time  step. 

For  a  time  step  i  the  loss  of  stress  is  expressed  as 


where 

CpsQ  =  thermal  strain  of  steel,  from  Eq.  (3.19) 
€<,e  =  thermal  strain  of  concrete,  from  Eq.  (3.10) 
6311  =  shrinkage  strain,  from  Eq.  (3.11) 
Ccr  =  creep  strain  of  steel,  from  Eq.  (3.20) 


(5.7 


.-.-^ 
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Epg  =  modulus  of  steel  at  the  beginning  of  time  step, 
from  Eq.  (3.4) 
The  effective  prestress  force  for  time  step  i  and  iteration 
j  is  confuted  by  following  expressions:  > 


f J  .   =  f^       -a 


(5.8) 


^1  ^ps,  1  "■ 


PS 


where 

Ap,  =  section  area  of  prestress  steel, 
o  =  cumulative  stress  loss,  from  Eq.  (5.7) 

5.4.  Evaluation  of  Action  Load  Vector 

The  action  load  vector  is  defined  as  the  effective    :  i 
loading  of  the  system.   It  is  derived  from  externally 
superinposed  node  and  element  loads,  equivalent  element  node 
loads  due  to  prestress  action  and  induction  forces  discussed 
in  Section  5.2.   The  transformation  of  the  in5>osed  element 
load  into  equivalent  node  loads  is  shown  in  Fig.  5.2.  ■ 

Equivalent  nodal  load  vector  is  then  expressed  as 

rI  =  {Rl,    R2,    R3,    R4,    R5 ,    R6} 

(5.9)   • 

=  ^{0,  6,  L,    0,  6,  -L)  -^ 

where  by  is  the  uniform  superimposed  load  acting  in  the 
positive  y-direction.   The  action  load  vector  is  the 
summation  of  results  of  Eqs .  (4.48),  (5.1),  (5.9) 
and  external  node  loads.   Thus  for  the  time  step  i 


:j 
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/f^|||/|v||)^;^/^/|^/^/|v,|i|;^j^/|>|j^/^ 


>,  R  4 


Fig.  5,2.  Equivalent  nodal  forces 


»  r 


and  iteration  j  the  action  load  vector   R  is  given  by 


^a,i*l    ~   ^n    '^   ^P,i*l    '*'   ^g,  i*l 


:5.io; 


where  vector  R^  is  induction  load  vector  given  by  Eq. (5.1) 
and  Rp  is  the  equivalent  prestress  load  vector  given  by 
Eq. (4.48)  which  have  been  adjusted  by  Eq. (5.8) .    R„  is 
constant  and  includes  elements  of  Eq. (5.9)  and  external 
point  loads.   Vector  R^^  is  computed  for  each  iteration  due 
to  the  changing  vectors  R^  and  Rp 

5.5.  Solution  of  Displacement  Response 

The  equilibrium  equations  to  be  solved  for  structural 
displacement  responses  are  Eqs.  (5.1)  and  (5.2).   Action 
load  vector  R  is  con5>uted  by  Eq.  (5.10).   Fig.  1.2.  shows 
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that  the  internal  load  function  F  attains  a  different  curve 
for  each  time  step  as  the  system  attains  a  new  deformed 
equilibrium  state  r.   This  characteristic  is  caused  by 
changes  in  prestress  load  function  P  and  stiffness  function 
S  as  depicted  in  Fig.  1.3.   The  shifting  behavior  of  the 
function  F  cuirve  presents  corr^jlications  in  the  accurate 
determination  of  the  new  value  of  F  at  the  beginning  of  each 
time  step. 

To  avoid  the  problems  of  tracing  a  new  curve  of  the 
function  F  in  each  time  step,  a  method  of  unbalanced  load 
adjustment  is  utilized.   The  strategy  in  this  approach  is  to 
determine  an  equivalent  incremental  load  function  which  will 
cause  increment  in  the  structure  deformation  from  the 
equilibrium  deformation  attained  in  the  previous  time  step. 
The  confutation  techniques  for  the  first  time  step  are 
therefore  different  from  the  other  time  steps.   The  solution 
process  in  the  first  time  step  is  shown  in  Fig.  5.3.         ^ 

The  initial  operations  involved  in  solving  Eqs.  (5.1) 
and  (5.2)  in  the  first  cycle  of  the  first  time  step  are 

(a)    r°  =  0 

ib)  K°  =  f{r°) 

(5.11) 
(c)    R^  =  f{L,P{r°) 

id)  K°Arl-R°  ;- 

where  letter  f  indicates  functional  dependence.   Operations 
in  subsequent  iterations  are  expressed  as  follows : 
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(a) 
ib) 
(c) 
(d) 
(e) 
if) 


K^  =  f{z^') 
AFi''  =  K^  hz^ 
R^  =  f{L,  P{z^^)) 
Ui'   =  Ri'   -   F/    ; 
K{Az^^  =  ui 


(5.12) 


R  ,   F 


Ar  1 
t  =  0 


t  =  0 


t  =0 


^  =  0     \-_, 


t  =  0 


s':^  r 


Fig.  5.3.  Solution  of  equilibrium  equations:  1st  time  step 


where  U  is  the  unbalanced  load  vector  to  be  determined. 
The  solution  process  for  other  time  steps  involves 
determination  of  incremental  action  load,  and  the 
incremental   internal  load  vector.  The  process  is 
presented  graphically  in  Fig .  5.4. 
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The  incremental  loads  to  be  confuted  at  the  beginning 
of  each  time  step  are 

(a)  induction  forces  due  to  thermal  and  shrinkage 
strain,  material  degradation  and  shift  in  stress-strain 
relation  law  and 

(b)  changes  in  the  prestress  action. 


I  +  1 


Fig.  5.4.  Solution  process  at  subsequent  time  steps 


The  operations  required  in  this  stage  to  compute  the 
incremental  loads  are  given  by  the  following  expressions: 


^■:>*-'.  ;-■ 
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^'-  (a)     rl^-  ri 

Ity  (^)  KU  =  f[j:l^) 

(c)  Aaf.i  =  o^(ri,e^)  -  a^,^[T^,^,e^) 

id)  G,,i  =  f(Aof.i,AoJ.i,Ao?;?)  (5.13) 

(e)  /^Rl,  =RljL,Pl,)    -R{{L,Pi) 

ig)  Kl^Arl^  =  ARl, 

where 

U  =  unbalanced  load  function  carried  forward  from  the 
previous  time  step, 
i'  S  =  indicates  induced  stress  due  to  shift  in  stress, 

T  =  indicate  induced  stress  due  to  thermal  effects 
(shrinkage  and  expansion) 
>  D  =  indicate  induced  stress  due  to  material 

■  '"  degradation  (crushing  or  cracking) 

Parameters  to  be  coit^juted  in  the  subsequent  iteration 
include 

(a)  current  deformation  of  the  system, 

(b)  update  of  stiffness  matrix, 

(b)  induction  force  due  to  material  degradation, 

(c)  update  of  prestress  action  and 

(d)  incremental  internal  load  function. 
Mathematically,  these  operations  are  described  in 

Eq.  (5.14).   The  total  incremental  displacement  within  a 
time  step  is  given  by  Eq.  (5.15) .   For  both  equations  j 
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(a)  r/.i  =  r/^t  +  5^/li 

(c)  6F/I1  =  ii:/li6r/li 

(d)  AF/li  =  Af/Ti^  +  6F/:i 

(e)  :»/,!  =  f  (L,  P/li) 

(5.14) 

if)  ^u  =  Ri,  -  Ri:i 

(ST)  bGf:i  =  f{oi,.,  (r/lj  ,  n^i.i  (r,.j ) 

{h)  Aie/li  =  Ai?/.;^  +  6rU  +  5c?/:i 

(i)  U/li  =  Ai?/li  -  AF/li 


Ari.,  =f5ri.,  ^^'^^^ 


represents  the  number  of  iterations  and   6ri+i  is  the 
displacement  increment  for  each  iteration.   Total  structure 
displacement  response  vector  at  the  end  of  iteration  j  for 
the  time  step  i+1  is  determined  by 

r1     =  r°-      +  Ar=!  (5.16) 

At  the  end  of  each  iteration  the  system  is  tested  for 
load  convergence.   The  acceptable  load  tolerance  level 
Ck  is  determined  by  the  following  expression: 


uL 


^{ArL   +  AF/li) 


<e  (5.17; 


83 


5.7.  Local  and  Global  Values 


The  con^jutations  above  involve  transformation  of 
matrices  and  vectors  from  local  element  to  global  structural 
coordinates,  and  vice  versa.   The  transformation  matrix   a  is 
expressed  as 


a  = 


c  s  0 

-SCO 

0  0  1 


0 

' 

c 

s 

0 

-s 

c 

0 

0 

0 

1. 

(5.18) 


where 

c  =  cos  9 

s  =  sin  9 

9  =  angle  between  global  and  local  x-axes 
Thus  the  element  stiffness  matrix  K^,  given  by  Eq.(4.35),  is 
transformed  to  structure  matrix  K,  by 


IC,  =  a^K.a 


(5.19) 


■'-1 


The  element  load  vector  (  R^)  due  to  superinposed 
uniform  load  given  by  Eq.  (5.9) ,  and  the  equivalent 
prestress  nodal  action  given  by  Eq.  (4.48),  is  transformed 
to  global  load  vector  R,  by  the  following  expression: 


^^ 


a^R^ 


(5.20; 


i 


-«  : 
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The  global  nodal  displacements  given  by  Eq.  (5.16)  are 
transformed  to  element  nodal  displacements  by 


5.8.  Confutation  of  Moment  Capacity  Response 

The  basic  moment  capacity  equation  for  flexural  members 
is  given  by  Eq.  (1.1) .   The  moment  equation  is  rewritten  as 


R  = 


^s   = 


(5.22) 


M^   =  t^l^sP'M^'  -   0.5i?r,FpLe)] 


1=1 


^ 


;^ 


where  '   *^ 

r^  =  terrperature  coefficient  of  coitpressive  strength  of 

concrete^  '  ><■■>    *•;> 

r^  =  coefficient  for  steel  stress  at  ultimate  load 
n  =  number  of  steel  segments  contained  in  a  beam 
element 
Within  the  limits  of  the  heat  transmission  end  point  f 

criteria  the  average  loss  of  compressive  strength  of  r  i 

f?  '■! 

concrete  in  the  corrpression  zone  is  2%  [6]  .   Therefore  the 
coefficient  r  ^  equals  0.98.   At  fire  teitperatures  the 
coefficient  r ,  is  taken  as  equal  to  0.98  [57]. 


CHAPTER  6 
COMPUTER  IMPLEMENTATION 


6.1.  Introduction 


The  numerical  procedures  presented  in  Chapters  2  to  5 
have  been  incorporated  in  the  conputer  program  PRECET        « 
(PREstressed  Concrete  subjected  to  Elevated  Tenyeratures) . 

This  chapter  describes  the  structure  of  the  program.  The 

■  ■■'■..     ■  ■  ^ 

instructions  for  preparation  of  input  data  file  are  given  in 
Appendix  A.  - 

6.2.  Main  Program  PRECET  -f  -      „  -"a 

The  program  PRECET  consists  of  all  the  routines  for 
reading  and  initial  processing  of  input  data.  As  the  primary 
executor  PRECET  transfers  initial  controls  to  data  modules 
and  then  to  executors  for  teirperature  analysis,  TEDIAN  and        :  v' 
for  response  analysis,  STREAN.   The  hierarchical  structure 
of  the  program  PRECET  is  presented  in  Fig.  6.1. 

6.3.  Program  TEDIAN 

TEDIAN,  TEmperature  Distribution  ANalysis,  is  the 
executor  for  temperature  distribution.  It  incorporates 
procedures  discussed  in  Chapter  2. 
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In  each  time  step  TEDIAN  is  called  first  to  determine 
temperature  distribution  across  the  beam  cross-section. 
Teirperature  values  for  the  subelement  are  saved  for  use  in  ' 
confuting  temperature-dependent  material  properties.  The 
flow  chart  of  the  program  TEDIAN  is  shown  in  Fig.  6.2. 

6.4.  Program  STREAN 

The  program  STREAN,  STructure  REsponse  ANalysis  is  the 
main  executor  of  response  analysis.   It  incorporates 
procedures  described  in  Chapters  4  and  5. 

STREAN  is  activated  at  each  time  step  after  ten^ierature 
analysis.   It  calls  the  material  property  routines  to     .  iJ 
compute  relevant  properties  of  each  subelement  based  on  its 
ten5>erature .  -       v; 

The  flow  chart  of  the  execution  logic  of  program  STREAN 
is  presented  graphically  in  Fig.  6.3,      *    ,  . 

6.5.  Program  MAPRAN 

MAPRAN,  MAterial  PRoperty  ANalysis,  consists  of  several 
independent  modules  each  of  which  confutes  a  specific 
property.  The  formulation  of  the  routines  is  based  on  the 
material  properties  and  procedures  presented  in  Chapter  3 . 

Thermal  property  routines  are  called  from  TEDIAN  while 
deformation  and  mechanical  property  modules  are  called  from 
program  STREAN. 
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Fig.  6.1.  Flow  Chart  of  Program  PRECET 
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Fig.  6.2.  Flow  Chart  for  Temperature  Distribution 
Analysis   Program  (TEDIAN) 
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Fig.  6.2.  Flow  Chart  for  Temperature  Distribution 
Analysis   Program  (TEDIAN)  (cont'd) . 
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Fig.  6.3.  Flow  Chart  for  Structure  Response 
Program  (STREAN) 
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CHAPTER  7 
NUMERICAL  ANALYSIS 


7.1.  Introduction 


This  chapter  presents  the  results  of  the  application  of 
program  PRECET  in  predicting  the  structural  response  of 
prestressed  concrete  members  subjected  to  elevated 
tender a ture s .   The  structural  configurations  analyzed  have 
been  subjected  to  standard  fire  tests.   The  results  of  the 
numerical  model  are  conpared  with  those  obtained  from  the 
fire  tests . 

Two  sets  of  problems  analyzed  are  prestressed  concrete 
beams  and  prestressed  concrete  slabs.   Results  presented 
include  teitperature  distribution,  mid-span  deflection  and 
%4l  i,    moment  capacity  ratio  for  beams,  mid-span  deflection  and 
moment  capacity  ratio  for  slabs. 


7.2.  Analysis  of  Prestressed  Concrete  Beams 

7.2.1.  Description  of  Test  Specimen 

The  beam  configurations  used  in  the  analysis  are  based 
on  the  fire  tests  conducted  on  eighteen  prestressed  concrete 
beams  at  the  Portland  Cement  Association's  Beam  Furnace  [2] . 
All  the  beeuns  had  identical  cross-section  and  of  the  same  ' fi': 

\ 
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span  of  20  ft.   Types  of  concrete  used  were  normal  weight 
and  lightweight  aggregates.   Three  cover  thicknesses  of  1, 
2,  and  3  in.  were  investigated  for  each  concrete  type. 
-^    The  specimens  were  divided  in  six  groups  each  of  three 
identical  beams.   The  first  three  groups  were  of  beams  made 
of  normal  weight,  silicious  aggregate  concrete  and  the  rest 
were  made  of  lightweight  expanded  shale  concrete.  The  three 
groups  in  each  concrete  type  were  categorized  in  terms  of 
the  concrete  covers  described  above . 

7.2.2.  Beam  Testing 

Fig.  7.1  (a)  shows  a  schematic  of  the  test  setup  for 
the  beam  fire  test.   During  the  test  steel  teitperatures  at 
selected  locations  in  concrete  were  measured  at  five 
sections,  A,  B,  C,  D  and  E  along  the  beam  length.   The  beams 
were  loaded  with  a  design  service  load  distributed  along  the 
beam  length  at  five  locations  as  shown  in  Fig.  7.1  (a) . 

Reported  test  results  for  each  group  have  been  averaged 
from  results  of  the  three  beams  in  each  group.   In  addition 
teirperatures  averages  were  determined  from  readings  along 
the  beam  length.   Tenperature  readings  from  the  extreme 
locations,  A  and  E,  were  lower  than  at  sections   B,  C  and  D. 
For  this  reason  only  ten^ierature  values  from  sections  B,  C 
and  D  were  used  in  confuting  the  averages . 
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Fig.    7.1.    Design  details  of  PC  beam  specimens 
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7.2.3.  Design  Details 

The  details  of  section  design  and  prestress  steel 
arrangements  for  the  three  beam  categories  are  shown  in 
Fig.  7.1  (c) .   Design  specifications  for  materials  and 
loadings  are  described  below. 
7.2.3.1.  Materials 


kC  > 


Prestressing  steel  specifications  were 

(a)  Ultimate  Strength,  f  p„  =  250,000  psi. 

(b)  Initial  Prestress,  f  ^i  =  0.7fp„  =  175,000  psi. 

(c)  Size  of  Prestress  Steel  Strand,   4),  =  3/8  in. 

(d)  Size  of  Prestress  Steel  Wire,  ^   =  0.195  in. 

(e)  Initial  Force  on  each  Strand  =  14,000  lb. 

(f)  Initial  Force  on  each  Wire  =  5,250  lb. 

(g)  Estimated  Loss  of  initial  Prestress  =  18%. 
Concrete  data  specifications  were 

(a)  Conpressive  Strength,  f '  „  =  5000  psi, 
at  28  days. 

(b)  Minimiim  Concrete  Strength,  f '  ^^  =  4000  psi,  at 
the  time  of  transfer  of  steel  stress  to 
concrete . 

(c)  Unit  Weight  of  Concrete  is  145  pcf  for  normal 
weight  and  97  pcf  (2/3  of  normal)  for 
lightweight. 
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7.2.3.2.  SuperinposGd  loads 

The  superinposed  load  at  each  loading  point  for  normal 
weight  beams  and  each  cover  category  are 

(a)  Beam  #  NWBEAMl:  P  =  28,000  lb. 

(b)  Beam  #  NWBEAM2 :  P  =  21,650  lb. 

(c)  Beam  #  NWBEAM3 :  P  =  12,500  lb. 
Loads  on  lightweight  beams  are 

(a)  Beam  #  LWBEAMl :  P  =  28,000  lb. 

(b)  Beam  #  LWBEAM2:  P  =  21,650  lb. 

(c)  Beam  #  LWBEAM3:  P  =  12,500  lb.  ■] 

7.2.3.3.  Other  design  parameters  '| 
The  minimum  age  of  concrete  during  fire  tests  was  one 

year.  During  that  period  all  the  beams  were  exposed  to  an 

environment  of  35%  RH  at  73  °F.  All  beams  were  subjected  to 

fire  test  while  supporting  their  design  service  loads  over  i 

a  span  of  20  ft.   The  limit  of  mid  span  deflection  was 

selected  at  15  in.   This  was  estimated  to  occur  only  a  few 

minutes  before  actual  structural  failure. 

7.2.4.  Numerical  Model 


Both  normal  and  lightweight  beam  specimens  were 
analyzed  for  terrperature  distribution,  displacement  response 
and  moment  capacity.   The  model  for  thermal  analysis 
consisted  of  a  mesh  of  116  nodal  points  and  176  triangular 
elements  over  the  half  width  of  the  cross  section. 
Structure  analysis  problem  was  modeled  using  a  system  of  11 
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beam  elements  and  12  nodes  over  the  half  beam  span.   The 
finite  element  idealization  of  the  structure  model  is  shown 
in  Fig.  7.1  (b) . 

7.2.5.  Numerical  and  Test  Result 


Numerical  results  for  temperature  distribution  and 
midspan  deflections  of  all  specimen  groups  have  been  plotted 
together  with  test  results  for  conparison.   No  test  results 
were  available  for  moment  capacity  ratio  therefore  only 
results  from  numerical  analysis  are  plotted. 
7.2.5.1.  Temperature  distribution  results 

Results  for  tenperature  distribution  for  certain  tendon 
locations  are  plotted  in  Figs.  7.2  to  7.4  for  normal  weight 
concrete  beams  and  Figs.  7.5  to  7.7  for  lightweight  concrete 
beams.   Numerical  results  show  relatively  good  agreement 
with  the  test  results  for  most  of  the  exposure  periods. 

Fig,  7.8  shows  average  temperature  distribution  of  the 
beam  sides  which  are  exposed  to  fire.   Tenperature  values 
for  lightweight  concrete  tend  to  be  slightly  higher  than 
those  for  normal  weight  concrete.   This  may  be  due  to  the 
relative  poor  ten5>erature  distribution  properties  of 
lightweight  concrete. 

Fig.  7.9  presents  the  average  coefficients  of  radiation 
heat  transfer  for  both  normal  weight  and  lightweight 
concretes.   The  computations  of  the  coefficients  are  based 
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Fig.  7.2.  Tenperature  distribution,  normal  weight 
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on  Eq.  (2.5).   Data  presented  in  Figs.  7.8  and  7.9  can  be 
useful  for  thermal  analysis  in  general  programs  like 
ANSYS  [58]  which  require  such  data  as  input. 
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Fig.  7.8.  Average  of  tenperatures  of  beam 
sides  exposed  to  fire 
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coefficients  for  beam  sides  exposed 
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7.2.5.2.  Mid  span  deflection  of  PC  beams 

Mid  span  deflection  responses  are  plotted  in  Figs.  7.10 
to  7.12  for  normal  weight  concrete  and  in  Figs.  7.13  to  7.15 
for  lightweight  concrete.   The  figures  show  that  the 
numerical  model  is  stiffer  than  the  test  model  for  most  of 
the  exposure  times.   The  reason  for  this  may  be  partly 
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Fig.  7.10.  Midspan  deflection,  normal  weight 
concrete  beam,  1  in.  cover 


attributed  to  the  assunptions  made  on  bond.   In  this  study  a 
perfect  bond  between  concrete  and  prestress  steel  is  assumed 
for  all  the  tenperature  ranges.   This  assumption  results  in 
full  transfer  of  prestress  action  from  steel  to  concrete. 
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7.11.  Mid  span  deflection,  normal  weight 
concrete  beam,  2  in.  cover 
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Fig .  7.12.  Mid  span  deflection ,  normal  weight 
concrete  beam,  3  in.  cover 


104 


ij-  \ 


n:;\ 


15 
14 
13 
12 
1  1 
10 

9  -| 

8 

7 

6 

5 

4 

3 

2 

1 

0 


1 :  Tesf  Results 

2:  Numerical  Results 


0  10         20         30         40         50         60         70 

Exposure   Time,     min 


80         90        100       110 
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concrete  beam,  1  in.  cover 


Fig.  7.14.  Mid  span  deflection,  lightweight 
concrete  beam,  2  in.  cover 
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Fig,  7.15.  Mid  span  deflection,  lightweight 
concrete  beam,  3  in.  cover 


concrete.   However,  according  to  tests  by  Fagundo  [19]  high 
terrperatures  cause  bond  slip  between  concrete  and  prestress 
steel.   Thus  in  the  real  structure  efficiency  of  transfer  of 
prestress  action  decreases  at  high  ten^serature .   This  may  be 
manifested  by  higher  deflection  results  for  the  test  model. 
On  the  other  hand,  numerical  results  tend  to  be  higher 
than  test  results  at  higher  tenperatures  close  to  the  end 
point.   The  reason  for  this  may  be  due  to  the  fact  that  at 
higher  teitperatures  significant  portions  of  the  section  have 
degraded.   The  algorithm  of  this  study  treats  all  the 
recently  degraded  elements  as  having  zero  stiffness  and 
carrying  no  stress.   In  reality  these  degraded  elements 
retain  some  diminishing  stiffness  and  stresses .   This 
discrepancy  results  in  a  higher  load  vector  and  lower 
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stiffness  values  for  the  numerical  model  which  is  then 
manifested  by  higher  deflections  for  the  numerical  model. 
7.2.5.3.  Moment  capacity  ratio 

Moment  capacity  ratio  is  defined  here  as  the  ratio  of 
the  service  load  moment  to  the  residual  teirperature- 
dependent  moment  capacity  of  a  member  considered.   The 
history  of  the  moment  capacity  ratio  for  normal  weight 
concrete  beams  is  presented  in  Fig.  7.16,  and  in  Fig.  7.17 
for  lightweight  concrete  beams . 
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Fig.  7.16,  Moment  capacity  ratios  for 

normal  weight  concrete  beams 


The  end  point  values  from  the  niimerical  models  of  this 
study  are  conpared  with  values  con5>uted  by  the  method 
described  by  Eq.  (1.1)  and  reported  in  Ref.  [6].   The 
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Fig.  7.17.  Moment  capacity  ratios  for 
lightweight  concrete  beams 
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Table  7.1.  Moment  Capacities  for  PC  Beams 


1 

2 

3 

4 

5 

Beam 
Type 

Concrete 
Type 

Concrete 
Cover, in. 

Moment  Capacity, in. ksi 

Eq. (1.1) 

Eq. (5.22) 

NWBEAMl 

NW 

1 

976 

1099 

NWBEAM2 

NW 

2 

778 

792 

NWBEAM3 

NW 

3 

492 

494 

LVTBEAMl 

LW 

1 

941 

1026 

LnBEAM2 

LW 

2 

744 

811 

LWBEAM3 

LW 

3 

458 

489 
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variables  used  in  Eq.  (1.1)  are  derived  from  measured  steel 
temperatures  at  the  end  of  test  time.   The  results  for  both 
numerical  procedures  are  shown  in  Table  7.1.   The  good 
agreement  of  the  results  shows  that  the  developed  numerical 
method  can  predict  well  the  fire  endurance  of  beams. 

7.3.  Prestressed  Concrete  Slabs 


Fire  endurance  tests  on  simply  supported  prestressed 
concrete  slabs  have  reported  by  Gustaferro  [21]  .   The  fire 
tests  were  conducted  on  eleven  pretensioned  concrete  slab 
specimens.   Six  of  the  specimens  were  made  of  normal  weight         j 
concrete,  and  the  rest  were  made  of  lightweight  concrete.  J; 

Due  to  the  qualitative  similarity  of  results  only 
conparison  results  of  normal  weight  concrete  slabs 
are  presented  in  this  section.  :i 


7.3.1.  Description  of  Test  Specimen 

All  the  slabs  used  in  this  series  of  tests  were  6.5  in. 
thick.   Three  thicknesses  of  concrete  cover  were  used,  that 
is,  1,  2,  and  3  in.   Three  of  the  normal  weight  concrete 
sleQss  were  pretensioned  with  five  7/16-in.  diameter 
seven-wire  strands.   The  cross  section  details  for  these 
sleQss  is  shown  in  Fig.  7.18  (b)  .   The  other  three  were 
pretensioned  with  fifteen  1/4-in.  diameter  seven-wire 
strands.   Fig.  7.19  (b)  shows  the  cross  section  details  of 
steel  arrangement  for  this  group. 
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c)  Finite  Element  Beam  Element  Model 


Fig.  7.18.  Design  details  for  slab  specimens, 
5-  7/16"  steel,  20  ft.  span 
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c)  Finite  Element  Beam  Element  Model 


Fig.  7.19.  Design  details  for  slab  specimen, 
15-  1/4"  steel,  12  ft.  span 
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7.3.2.  Material  Specifications 


The  measured  ultimate  tensile  strength  for  the  1/4-in. 
strands  was  275,000  psi  and  262,000  psi  for  the  7/16-in. 
strands .   All  the  strands  were  initially  pretensioned  to  a 
nominal  175,000  psi.   Partial  loss  of  prestress  is  assumed 
to  be  18%.   Details  of  concrete  strengths  and  loadings  are 
presented  in  Teible  7.2. 
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Table  7.2.  Specimen  strength  and  loading  details. 

(Source:  Ref.  [21]) 


Slab 
Type 

Concrete 
Cover ,  in . 

Concrete 
Strength 
(psi) 

Number  and 
Size  of 
Strands 

Applied 
Load,  P 
(lb) 

NWSLABl 

1 

5400 

5-7/16  in. 

1297 

NWSLAB2 

2 

5010 

5-7/16  in. 

977 

NWSLAB3 

3 

6300 

5-7/16  in. 

643 

NWSLAB7 

1 

5180 

15-1/4  in. 

2569 

NWSLAB8 

2 

7740 

15-1/4  in. 

2035 

NWSLAB9 

3 

5430 

15-1/4  in. 

1365 

7.3.3.  Testing  Procedures 


All  specimens  were  singly  supported  on  rocker-rollers, 
and  fire  tested  while  supporting  their  designed  service 
loads.   The  loading  was  applied  through  five  equally- spaced 
hydraulic  rams,  as  shown  in  Figs.  7.18  (a)  and  7.19  (a). 
Specimens  pretensioned  with  7/16-in.  strands  were  tested  on 
a  20-ft.  span.   The  other  specimens  were  tested  on  12-ft. 
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span.   The  relative  humidity  of  the  specimens  at  the 
beginning  of  the  fire  test  was  70%.   Tests  were  conducted 
until  the  structural  end  point  was  imminent. 

The  fire  exposures  were  conducted  in  accordance  with 
ASTM  E-119.   Ten^erature  of  steel  strands  was  measured  by 
attaching  thermocouples  to  several  of  the  strands  near 
midspan  and  at  a  distance  of  5-ft.  from  midspan. 

7.3.4.  Finite  Element  Model 


The  finite  element  model  for  thermal  analysis  consisted 
of  224  triangular  elements,  and  135  node  points  over  one 
half  of  the  slab  width.   For  structure  response  analysis  the 
half  span  length  of  the  specimens  was  modeled  int  10  beam 
elements.   The  idealization  is  shown  in  Fig.  7.18  (c)  for 
the  20-ft.  span  and  in  Fig.  7.19  (c)  for  the  12-ft.  span. 

7.3.5.  Conyarison  of  Numerical  and  Test  Results 

7.3.5.2.  Temperature  distribution 

Fig.  7.20  shows  the  average  tenperature  of  the  face  of 
the  slab  exposed  to  fire  as  function  of  time.   Also 
presented  is  the  average  and  the  maximum  temperature 
histories  of  the  unexposed  surface.   The  curve  for  the 
coefficient  of  radiation  heat  transfer  for  normal  weight 
concrete  sledss  is  presented  in  Fig.  7.21.   Conparison  of 
numerical  and  test  results  for  teitperature  distribution 
across  slab  section  are  not  presented  because  the 
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corresponding  test  results  were  not  availedsle.  As  for  the 
beams  the  exposed  surface  teirperatures  given  in  Fig.  7.20 
and  the  radiation  heat  transfer  coefficients  given  in 
Fig.  7.21  can  be  used  in  general  thermal  analysis  programs 
as  input  data  for  confuting  teirperature  distribution  in 
slabs  subjected  to  fire  tests. 
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Fig.  7.20.  Exposed  and  unexposed  surface  temperatures, 
normal  weight  concrete  slabs 
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Fig.  7.21.  Coefficient  of  radiation  heat  transfer 
for  normal  weight  slab  exposed  surfaces 
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The  curve  for  unexposed  surface  temperature  and  the 
maximum  tenperature  curve  determine  the  heat  transmission 
end  point  for  slabs . 
7.3.5.2.  Mid-span  deflection  of  slabs 

The  coirparison  of  test  and  numerical  results  of  mid 
span  deflections  for  normal  weight  concrete  slabs  is 
illustrated  in  Figs.  7.22,  7.23  and  7.24  for  20-ft  span,  and 
in  Figs.  7.25,  7.26  and  7.27  for  the  12-ft  span.  The 
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Fig.  7.22.  Mid  span  deflection,  normal  weight  slab, 
5-  7/16"  steel,  1  in.  cover,  20  ft.  span 


results  in  these  figures  shows  only  the  increase  of 
deflection  due  to  effects  of  the  increase  in  temperature. 
The  initial  deflection  due  to  loading  and  prestress  action 
has  been  deducted  from  the  computed  deflection  results. 
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As  was  the  case  with  the  beams  numerical  results  are 
relatively  lower  than  the  corresponding  test  results  for 
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Fig.  7.23.  Mid  span  deflection  for  normal  weight 
slab,  5-  7/16"  steel,  2  in.  cover, 
20  ft.  span 
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7.24.  Mid  span  deflection,  normal  weight  slab, 
5-  7/16"  steel,  3  in.  cover,  20  ft.  span 
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Fig.  7.25.  Mid  span  deflection,  normal  weight  slab, 
15-  1/4"  steel,  1  in.  cover,  12  ft.  span 
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1 :  Test  Results 


2:  Nunnerical  Results 
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Fig.  7.26.  Mid  span  deflection,  normal  weight  slab, 
15-  1/4"  steel,  2  in. cover,  12  ft.  span 
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Fig.  7.27.  Mid  span  deflection,  normal  weight  slab, 
15-  1/4"  steel,  3  in.  cover,  12  ft.  span 
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most  of  the  lower  tenperatures .   It  can  also  be  noticed  that 
the  discrepancy  between  niimerical  and  test  results  is 
greater  for  slabs  than  for  beams.    However,  the  arguments 
on  the  results  presented  for  beams  are  also  valid  for  slabs . 
7.3.5.3.  Moment  capacities  ratios  for  slabs 

The  variations  of  moment  capacity  ratios  as  related  to 
exposure  time  are  presented  in  Fig.  7.28  for  the  20-ft.  span 
slabs  and  in  Fig.  7.29  for  the  12-ft.  span  slabs.   The  end 
point  moment  capacities  generated  from  the  method  of  this 
study  and  those  of  the  method  of  Eq.  (1.1)  are  conf>ared  in 
Table  7.3.   The  closeness  of  results  verifies  the  validity 
of  the  niunerical  model  used  in  this  study. 
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Fig.  7.28.  Moment  capacity  ratios  for  normal 

weight  concrete  slabs  of  20  ft.  span 
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Fig.  7.29.  Moment  capacity  ratios  for  normal 
weight  slabs  of  12  ft.  span 
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Tedsle  7.3  End  point  moment  capacities  of  PC  slabs 
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Slab 
.-  Type 

Span 
ft. 

Concrete 
Cover, in. 

Moment  Capacity, in. ksi 

Eq.(l.l) 

Eq.  (5.22) 

NWSLABl 

20 

1 

336 

352 

NWSLAB2 

20 

2 

315 

336 

NWS LABS 

20 

3 

214 

227 

NWSLAB7 

12 

1 

283 

329 

NWSLAB8 

12 

2 

230 

250 

NWSLAB9 

12 

3 

178 

166 

7.4.  Evaluation  of  Influences  of  Material  Properties 

-*■■ 


Data  in  Chapter  3  demonstrates  the  effect  of  higher 
tenperatures  on  properties  of  concrete  and  prestressing 
steel.   These  effects  on  material  properties  have  different 
influences  on  the  teitperature  distribution  and  structural 
response  observed  above.   This  section  presents  the 
evaluations  of  influences  of  certain  material  properties  on 
the  computed  deflections,  moment  capacities  and  tenperature 
distributions . 

Nvimerical  results  based  on  the  tenperature  dependent 
properties  of  concrete  and  steel  are  considered  the  base  for 
evaluating  the  influences  of  the  effects  of  tenperature  on 
material  properties .   The  technique  used  in  determining  the 
influence  of  a  particular  parameter  is  to  treat  the 
parameter  as  independent  of  tenperature  during  the  analysis. 
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7.4.1.  Effect  of  Thermal  Properties  on  Tenyerature 

Distribution 


Thermal  conductivity,  thermal  capacity  and  their 
combinations  are  alternatively  considered  constant  in  the 
analysis  of  tenperature  distribution.   The  resulting  sets  of 
teirperature  distributions  are  plotted  in  Fig.  7.30  for 
normal  weight  concrete  and  in  Fig,  7.31  for  lightweight 
concrete.   The  curves  denote  the  tenperature  of  tendon 
location  marked  by  number  11  as  shown  in  both  figures. 
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Fig.    7.30.    Influences   of   thermal  properties   in  temperature 
distribution  of  normal   weight  beams 


For  normal  weight  concrete  the  maximum  tenperature 
value  is   increased  by  42%  when  both  conductivity  and  thermal 
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capacity  are  assumed  constant.   When  only  thermal 
conductivity  is  temperature  dependent  the  maximum  based 
tenperature  is  increased  by  27%,  while  the  increase  is  16% 
when  only  thermal  capacity  is  ten^erature  dependent.   This 
analysis  show  that  temperature  distribution  in  normal  weight 
concrete  is  more  sensitive  to  thermal  capacity  than  thermal 
conductivity. 
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Fig.  7.31.  Influences  of  thermal  properties  on  temperature 
distribution  for  lightweight  concrete  beams 


Thermal  analysis  for  lightweight  concrete  with  constant 
thermal  properties  results  in  24%  increase  in  the  maximum 
base  tenperature.   When  only  thermal  capacity  changes  with 
tenperature  the  base  tenperature  is  increased  by  20%,  and  by 
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6%  when  only  thermal  conductivity  is  tenperature  dependent. 
The  data  indicate  that  tenperature  distribution  for 
lightweight  concrete  is  highly  sensitive  to  the  thermal 
properties  of  the  concrete. 

7.4.2.  Effect  of  other  Material  Properties  on  Deflection 

In  this  analysis  the  thermal  properties  of  concrete  are 
tenperature  dependent.   Other  properties  are  grouped  and 
given  identification  symbols  as  follows:  modulus  properties 
(Ec,  Es) ,  strength  properties  (fc,  fpu) ,  deformation 
properties  which  include  thermal  and  shrinkage  strain  (Def.) 
and  relaxation  properties  of  steel  (Rel . ) .   The  parameters 
whose  influences  are  to  be  evaluated  are  assumed  to  be 
independent  of  tenperature.   The  results  will  reflect 
sensitivity  of  the  basic  deflection  on  the  evaluated 
parameters.   The  model  used  for  the  analysis  is  the  normal 
weight  beeun  with  2  in.  cover. 

The  deflection  curves  based  on  various  modes  of 
teitperature-dependence  material  properties  considerations 
are  presented  in  Fig.  7.32.   The  plotted  data  show  that  the 
influence  of  tenperature  on  individual  property  groups 
investigated  is  mostly  significant  only  at  higher 
tenperatures  near  the  structural  end  point. 

Evaluation  of  results  show  the  following.   When 
relaxation  property  of  steel  is  independent  of  tenperature 
the  time  taken  to  reach  the  15  in.  deflection  limit  is 


123 


15 
14 

13  H 
12 

^  1  1 

c 

■^  10  - 

c 

■B     9 
o 
®     8 


C 
0 
CL 
I/) 
I 

■q 

is 


7  H 

6 

5  - 

4 

3 

2 

1 

0 


props=f(t) 
def.=cons1 


Ec,  Eps=cons1 
Rel.=consl 


1 1 1 1 1 1 r 1 1- 


10   20   30   40   50 


60   70   80   90   100  110  120  130  140  150  160 

Exposure  Time.  min. 


Fig.  7.32.  Influence  of  non-thermal  material  properties 
on  deflection  for  normal  weight  beams 


increased  by  26%.   If  the  influence  of  teir?>erature  on 
deformation  properties  is  neglected  the  end  time  is 
increased  by  18%  while  the  increase  is  only  4%  for  constant 
strength  properties.   When  the  modulus  properties  are 
assumed  independent  of  temperature  the  resulting  deflection 
is  lower  than  that  of  other  groups  for  lower  teirperatures . 
This  is  probeibly  due  to  higher  modulus  values  which  increase 
structural  stiffness.   However  in  this  case,  instead  of 
increasing,  the  time  to  reach  end  point  deflection  is 
reduced  by  4% . 

This  analysis  shows  that  the  deflection  behavior  of 
prestressed  concrete  member  exposed  to  fire  is  mostly 
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sensitive  to  the  effects  of  teitperature  on  the  relaxation 
properties  of  the  prestressing  steel. 

7.4.3.  Effect  of  Thermal  Properties  on  Structural  Response 


This  section  demonstrates  how  thermal  properties  of 
^*'  concrete  affect  deflections  and  moment  capacities  of 

^  ^     prestressed  concrete  members  subjected  to  fire.   The 

analysis  is  carried  out  with  inproved  thermal  properties  of 
'       concrete.   In  this  analysis  thermal  properties  are  improved 
by  reducing  thermal  conductivity  values  and  by  increasing 
thermal  capacity  values.   Thus  5%  in5>rovement  means  the 
basic  temperature  dependent  thermal  conductivity  values  are 
reduced  by  5%  while  thermal  capacity  values  are  increased  by 
the  same  level.   All  other  material  properties  are 
considered  at  their  basic  functional  values. 

Four  levels  of  modifications  have  been  arbitrarily 
chosen  for  consideration:  5%,  10%,  15%  and  20%.   Deflections 
resulting  from  these  considerations  are  plotted  and 
presented  in  Fig.  7.33.   Corresponding  moment  capacity 
ratios  are  shown  in  Fig.  7.34. 

The  deflection  results  indicate  that  irt^jroving  thermal 
properties  by  only  5%  would  increase  the  end  point  time  by 
5%,  while  a  10%  iitprovement  result  to  a  19%  increase. 

Results  of  moment  capacity  show  that  a  5%  iit^irovement 
of  thermal  properties  increases  fire  endurance  by  7% .   An 
irrprovement  of  10%  increases  the  fire  endurance  by  16%. 
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Fig.  7.33.  Influences  of  thermal  properties  of  concrete 
on  deflections  of  normal  weight 
concrete  beams 
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Fig.  7.34.  Influences  of  thermal  properties  of  concrete 
on  moment  capacities  of  normal  weight 
concrete  beams 
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7.5.  Method  of  Estimating  Fire  Endurance  of  Beams 

7.5.1.  Effects  of  Exposure  Conditions 

^  «  •         Slabs  are  exposed  to  fire  on  the  bottom  surface  only. 
■^  '  '     Thus  the  teitperature  of  steel  in  the  slabs,  and  hence,  the 
sf\  %      fire  endurance  of  pres tressed  concrete  slabs,  can  be 

predicted  satisfactorily  using  the  numerical  procedures 
'    described  in  Sec.  1.6.3,  pg.  19  of  this  document  and  in 

Ref .  [3] .   The  accuracy  of  the  procedure  has  been  confirmed 
by  the  method  of  the  current  study  as  shown  in  Table  7.3. 

Beams,  however,  are  exposed  to  fire  on  three  surfaces. 
This  exposure  condition  makes  it  difficulty  to  estimate 
tenperature  of  prestressing  steels  using  the  tenperature 
curves  [3] .   Thus  the  numerical  procedure  described  in  [3] 
is  not  quite  accurate  for  predicting  fire  endurance  of 
prestressed  concrete  beams . 

This  section  presents  simple  design  charts  for 
predicting  fire  endurance  of  prestressed  concrete  beams. 
The  procedures  of  developing  the  charts  by  using  the  method 
of  this  study  is  also  presented. 

7.5.2.  Specimen  for  Analysis 

The  beam  specimen  used  in  the  analysis  is  shown  in 
Fig.  7.35.   The  number  and  arrangement  of  the  steel  tendons 
and  the  values  of  the  variables  depth  h,  width  b  and  tendon 
eccentricity  e  are  based  on  the  design  recommendations  given 
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in  page  2.24  of  Ref.  [57].   The  minimum  spacing  s  between 
steel  tendons  is  given  in  Sec.  7  of  Ref.  [32] .   This 
approach  is  utilized  in  order  to  produce  uniformity  in 
tendon  arrangement  for  all  the  cover  thicknesses  considered. 
Thus  the  variation  of  the  initial  moment  capacity  between 
beams  of  the  same  size  will  be  due  to  the  influence  of  cover 
thickness  on  the  number  of  tendons  that  can  be  accommodated 
in  the  beams. 
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Fig.  7.35.  A  typical  PC  beam  section 


7.5.3.  The  Analysis  Procedure  and  Results 


The  analysis  involved  three  cover  thickness  groups  of 
1,  2,  and  3  in.  in  normal  and  lightweight  aggregate 
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concretes.   The  design  variables  considered  in  each  cover 

group  are  beam  size,  tendon  strength  and  moment  intensity. 

The  moment  intensity  is  defined  here  as  the  ratio  of  service 

load  moment  M  to  the  moment  capacity  M  „  of  the  beeun  specimen 

at  normal  teitperature .   Beam  sizes  investigated  were  12  and 

16-in.  width,  and  24  and  32-in.  height.   The  tensile 

strength  of  tendons  were  250  and  270-ksi.   The  moment 

intensities  ranged  from  0.2  to  0.7. 

By  arranging  steel  tendons  in  layers  and  maintaining 

the  same  tendon  eccentricity  [57]  and  size  it  was  possible 

to  arrange  in  a  beam  the  same  number  of  tendons  for  both  1  ..-^ 

''i 
and  2-in.  covers  thicknesses.   However,  the  number  of 

tendons  were  less  for  3  in.  cover  thickness  due  to 

limitation  of  clearance  between  strands .   This  tendon 

arrangement  is  adopted  only  for  convenience  of  analysis.   It        ■{ 

does  not  have  a  significant  effect  on  the  resulting  curves 

-I 

since  fire  endurance  is  expressed  as  function  of  a 

''■i 

dimensionless  parameter,  that  is,  moment  intensity. 

For  each  cover  group  the  beam  specimen  involving  the 
design  variables  described  above  are  analyzed  for  fire  , 

endurance.   The  analysis  is  conducted  as  follows.   Consider, 
for  exeunple,  a  2  in.  cover  group.   For  this  group  a  beam  of 
size  12  by  24-in.,  tendon  strength  250  ksi  and  moment 
intensity  0.3  is  analyzed  to  determine  the  history  of  moment 

capacity  ratio.   Holding  the  variables  of  size  and  tendon 

■  i 
strength  constant  the  beam  is  analyzed  for  other  j 
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Fig.  7.36.  Moment  capacity  ratio  versus  time  for  a  normal 
weight  concrete  beam,  2  in.  cover  thickness 


moment  intensity  values .   The  results  of  all  the  analyses 
for  the  whole  range  of  moment  intensity  are  then  plotted 
together  as  shown  in  Fig.  7.36. 

The  fire  endurance  period  for  each  moment  intensity 
value  investigated  is  then  determined  as  the  time  when  the 
service  moment  M  equals  the  temperature-dependent  moment 
capacity  M^g.   These  values  of  fire  endurance  versus  moment 
intensity  form  a  set  of  results  of  fire  endurance  for  a  12 
by  24-in  beam  prestressed  with  tendons  of  250  ksi  strength. 
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This  process  is  repeated  to  produce  other  sets  of  results  of 
fire  endurance  for  all  the  combination  of  beam  sizes  and 
tendon  strengths  within  the  2  in.  cover  group. 

All  the  sets  of  results  of  fire  endurance  within  the 
2  in.  cover  group  are  plotted  together.   The  reason  for  this 
approach  is  to  determine  the  range  of  scatter  of  fire 
endurance  data  as  influenced  by  beam  size  and  tendon 
strength.   This  step  is  illustrated  in  Fig.  7.37.   The 
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Fig.  7.37.  Fire  endurance  of  various  beam  sizes  and  steel 
strengths  for  a  2  in.  cover  thickness, 
normal  weight  concrete 


whole  analysis  process  is  also  performed  for  1  and  3-in. 
cover  groups  and  for  both  normal  and  lightweight  concretes 


'i  t  >^ 


131 
The  set  of  fire  endurance  results  plotted  in  Fig.  7.37 
are  very  close  together.   This  indicates  that  beam  sizes  and 
tendon  strengths  have  no  significant  effect  on  the  fire 
endurance  of  a  prestressed  concrete  beam  with  2  in.  cover 
thickness.   Similar  behavior  has  also  been  observed  for 
beams  of  1  and  3-in.  concrete  covers.   The  final  values  of  • 

fire  endurance  of  prestressed  concrete  beeuns  are  confuted  by         , 
averaging  all  the  sets  of  fire  endurance  results  within  a 
cover  group.   The  resulting  values  of  fire  endurance  for  all 
the  cover  groups  considered  are  plotted  together  to  form 
design  charts.  ,j 

The  developed  design  charts  for  estimating  fire 
endurance  of  prestressed  concrete  beams  are  shown  in 
Fig.  7.38  for  normal  weight  concretes  and  in  Fig.  7.39  for 
lightweight  concretes.   The  charts  involves  the  variables  ',. 

concrete  cover,  moment  intensity  and  fire  endurance.   When         ■  | 
any  two  of  these  parameters  are  known,  the  third  can  be  | 

determined. 

The  application  of  the  charts  in  determining  fire 
endurance  is  demonstrated  through  the  following  exeunple.   It        ^' 
required  to  determine  the  fire  endurance  of  a  12  by  24-in. 
prestressed  concrete  rectangular  beam  (12RB24)  prestressed 
with  10-1/2-in.  diameter  tendons.   The  beam  section  and  the 
design  data  are  given  in  Fig.  7.40.   The  basic  calculations 
needed  are  for  the  determination  of  the  load  intensity. 
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Fig.  7.38.  Fire  endurance  as  function  of  moment 
intensity,  for  normal  weight  PC  beams 
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Fig.  7.39.  Fire  endurance  as  function  of  moment 
intensity,  for  lightweight  PC  beams 
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Fig.  7.40  Section  details  for  a  12RB24  beam 


The  solution  is  as  follows: 

Ap,  =  10(0.153)  =  1.53  in  2 
y,  =  [5(2.5)+5(4.5)]/10  =  3.5  in. 
dp  =  24-3.5  =  20.5  in. 

Using  Eq.  (1.1) 

fp.  =  270[1-(0.5*1.53*270)/(12*20.5*5)  =  224.7  ksi . 
a  =  (1.53*224.7)7(0.85*5*12)  =  6.74 

Nominal  moment  strength  is  given  by 

M„  =  (1. 53*224. 7)*(20. 5-6. 74)  =  4730  kips-in. 

Applied  moment  is 

M  =  wlV8   =  2.2(30)  V8  =  247.5   kips-ft. 
=  2970   kips-in. 

Load  intensity  is 

M/M„  =  2970/4730  =0.63 

With  M/M„=0.63  and  cover  c=2.5  in.  fire  endurance  of  the 

given  beam  is  determined  from  Fig.  7.38  as  2  hr  and  45  min. 
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CHAPTER  8 
SUMMARY  AND  CONCLUSIONS 


8.1.  Summary 


A  mimerical  procedure  and  computer  program  for  a 
nonlinear  analysis  of  the  structural  response  of  prestressed 
concrete  members  subjected  to  elevated  tenperatures  has  been 
presented.   The  parameters  defining  the  structural  response 
are  deflection  and  moment  capacity.   The  method  applies  the 
time- temperature  relationship  of  the  standard  fire  test 
(SFT)  to  simulate  building  fires.   The  numerical  model  is 
capable  of  determining  temperature  distributions  across 
prestressed  concrete  member  sections,  and  predicting 
deflections  and  moment  capacity  of  such  a  member  supporting 
its  designed  service  loads,   throughout  the  simulated  fire 
exposure  period. 

The  effects  of  elevated  teit^jeratures  on  the  material 
properties  of  concrete  and  prestressing  steel  are  included 
in  the  model.    The  material  properties  considered  are: 
thermal  conductivity  and  thermal  capacity  of  concrete, 
stress-strain  relationship  for  both  concrete  and 
prestressing  steel,  cortpressive  and  tensile  strength  of 
concrete,  tensile  strength  of  steel,  thermal  strain  of 
concrete  and  steel,  concrete  shrinkage,  and  relaxation  of 
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prestressing  steel.   The  behavior  of  these  properties  as 
related  to  teirperature  has  been  derived  from  data  published 
in  the  literature.   This  information  is  developed  into 
mathematical  relations  to  be  used  in  the  numerical  model. 

The  formulation  of  the  numerical  model  and  the  solution 
of  the  equilibrium  equations  is  based  on  the  displacement 
finite  element  method.   The  member  cross  section  is 
idealized  into  a  finite  number  of  triangular  elements  for 
temperature  distribution  analysis.   This  idealization  also 
is  used  to  account  for  the  variation  of  material  properties 
within  a  section.   The  properties  of  individual  elements  are 
numerically  integrated  across  the  section  in  order  to 
determine  the  section's  stiffness  properties.   The  member  is 
modeled  into  discrete  number  of  beeun  elements  for  the 
deflection  and  moment  capacity  analysis. 

The  exposure  time-dependent  analysis  is  performed  by 
dividing  the  time  scale  into  finite  time  step  intervals . 
For  each  time  step  equilibrium  equations  are  set  and  solved 
for  tenperature  distributions  and  structural  response.   The 
solution  of  the  equilibrium  equations  for  deflections  is 
iterative  and  includes  determination  of  the  incremental 
loads  induced  into  the  system  due  to  changes  in  temperature 
dependent  material  properties .   Incremental  deflections  for 
each  iteration  are  cumulatively  added  along  the  time  domain 
to  determine  total  deflection  for  each  time  interval. 
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The  con^uter  program  iirplementation  of  the  developed 
:^'i^'     numerical  procedures  is  used  to  analyzes  a  series  of 

prestressed  concrete  beams  and  slabs  which  have  been  fire 
tested  according  to  SFT  procedures .   The  tenperature 
distribution  and  structural  response  results  of  the 
numerical  model  are  conpared  graphically  with  those  obtained 
from  fire  tests  in  order  to  investigate  the  validity  and 
applicability  of  the  numerical  method.   Analyses  have  also 
been  performed  to  evaluate  the  degree  of  influences  of 
categories  of  temperature-dependent  properties  of  concrete 
and  steel  on  the  confuted  structural  response. 
The  numerical  model  has  also  been  used  to  develop  charts  for 
estimating  fire  endurance  of  prestressed  concrete  beams . 


8.2.  Conclusions 

Based  on  the  structural  response  analyses  performed, 
and  the  con5>arison  made  between  numerical  and  test  results, 
the  following  conclusions  are  reached  regarding  the 
numerical  procedures  developed  together  with  the 
corresponding  conputer  implementations: 

1.   The  numerical  model  is  capable  of  predicting 
deflection  response  of  prestressed  concrete  members 
subjected  to  elevated  teirperatures .   In  accordance  with 
results  of  fire  tests  utilized  in  the  evaluations,  the  model 
is  more  accurate  for  predicting  deflections  of  besuns 
than  of  sleU3s. 
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2 .  The  developed  numerical  method  is  also  capable  of 
predicting  the  history  of  the  moment  capacity  of  a 
prestressed  concrete  member  throughout  the  fire  exposure 
period.   The  estimation  of  the  structural  end  point  based  on 
moment  capacity  ratio  agrees  well  with  results  from  other 
previously  developed  numerical  approaches . 

3 .  The  computer  implementation  of  the  numerical 
procedures  developed  in  this  study  provides  means  of 
evaluating  the  sensitivity  of  the  coii5)uted  structural 
response  to  influences  of  tenperature  on  material 
properties . 

Evaluation  of  the  results  obtained  from  using  the 
current  numerical  procedures  to  determine  degrees  of 
influences  of  combinations  of  teirperature  dependent  material 
properties  leads  to  the  following  conclusions : 

1.  For  normal  weight  concretes,  tenperature 
distribution  is  more  sensitive  to  the  influence  of 
tert5>erature  on  thermal  capacity  than  on  thermal 
conductivity.   However,  tenperature  distribution  in 
lightweight  concretes  is  highly  sensitive  to  the  influence 
of  teirperature  on  thermal  conductivity. 

2 .  The  influence  of  ten?>erature  on  the  relaxation 
properties  of  prestressing  steel  have  the  highest  effect  on 
deflection,  followed  by  the  combined  effect  teit5>erature  on 
shrinkage  properties  of  concrete,  and  thermal  expansion 
properties  of  concrete  and  prestress  steel. 
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3 .   Slight  inprovements  in  thermal  qualities  of 
concrete  can  result  in  significant  inprovement  in  the 
performance  of  prestressed  concrete  members  at  high 
temperatures . 

8.3.  Recommendations 


In  accordance  with  the  conclusions  reached  above  and 
the  results  presented,  the  following  recommendations  are 
given : 

1 .  For  enhancing  the  performance  of  prestressed 
concrete  members  in  fire  environment  much  effort  should  be 
devoted  in  in^jroving  the  thermal  properties  of  concrete. 

2 .  More  studies  to  be  conducted  in  the  area  of  bond 
slip  and  its  effects  on  the  performance  of  prestressed 
concrete  members  subjected  to  high  tenperatures . 
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APPENDIX  A 
INSTRUCTIONS  FOR  THE  PREPARATION  OF  THE  INPUT 
DATA  FILE  FOR  THE  PRECET  PROGRAM 


INTRODUCTION 


This  document  presents  the  information  required  to 
create  an  input  file  for  the  program  PRECET.   The  document 
is  arranged  into  sections,  which  form  the  data  categories  of 
the  input  file.   The  input  data  consists  of  information  for 
both  the  thermal  and  the  structure  model . 

Due  to  symmetry,  only  half  of  the  cross-section  is 
considered  in  the  thermal  model.   The  section  is  modeled 
using  only  triangular  elements.   Each  element  must  have  one 
side  (base)  parallel  to  the  x-axis .   Elements  are  numbered 
in  counter-clock-wise  direction  in  the  i-j-k  sequence. 

The  structure  is  modeled  using  two  dimensional  beam 
elements.   The  uniform  distributed  loads  (element  loads)  are 
defined  positive  when  their  direction  is  towards  the 
positive  y-axis. 

In  entering  model  data,  two  back-slash  symbols,  '\\', 
have  to  be  used  to  indicate  continuation  of  data  covering 
more  than  one  line.  Total  number  of  characters  no 
to  exceed  160. 


SECTION  ONE:  PROBLEM  DEFINITIONS 
Line  #1 :  Problem  Heading 

PR0BLEMHEADIN6 
Line  #2:  Tenperature  dependence  control  data 

NRl  ,NR2  ,NR3  ,NR4  ,NR5  ,NR6  ,NR7  ,NR8  ,NR9  ,NR10  ,NR11  ,NR12 

where  the  entries  have  value  either  1  or  0  in  which 

NRI  =  0  means  property  not  tenperature  dependence, 
for  I  =  1  to  10;  thermal  analysis  only, 
for  I  =  11;  and  moment  capacity 
analysis  only,  for  I  =  12. 
=  1  means  property  is  tenperature  dependence, 
for  I  =  1  to  10;   thermal  and 
moment  capacity  analysis,  for  I  =  11; 
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and  coirplete  analysis,    for  I  =  12. 

The  properties  in  considerations  are 

NRl  =  thermal  conductivity  of  concrete 
NR2  =  thermal  capacity  of  concrete 

j^  <  .  ,   NR3  =  stress-strain  relation  for  concrete 

■^f^'  T.  ;  NR4  =  stress-strain  relation  for 

prestressing  steel 

--  ^  NR5  =  strength  of  concrete 

^<  -'-.i  (conpressive  and  tensile) 

NR6  =  tensile  strength  of  concrete 

'^»  NR7  =  thermal  strain  of  concrete 

^ • T '  NR8  =  thermal  strain  of  prestressing  steel 

NR9  =  shrinkage  strain 
NRIO  =  relaxation  of  prestressing  steel 

SECTION  TWO:  ITERATION  CONTROLS 

Line  #1 :  Convergence  control  header 

CONVERGENCE 

Line  #2 :  Convergence  criteria 

NRl ,NR2 ,RD3 ,RD4 ,RD5 ,RD6 ,RD7 

where , 

NRl  =  maximum  number  of  iterations  in 

thermal  analysis 
NR2  =  maximum  number  of  iterations  in 

response  analysis 
RD3  =  allowed  error  in  nodal  ten^eratures 
RD4  =  allowed  error  in  radiation  coefficient 
RD5  =  allowed  error  in  nodal  forces 
RD6  =  allowed  error  in  nodal  displacement 
RD7  =  critical  deflection 

SECTION  THREE:  MATERIAL  PROPERTIES 

Line  #1 :  Material  control  header 

MATERIALS 
Line  #2:  Property  data 

NRl , NR2 , RD3 , RD4 , RD5 , RD 6 
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where , 

NRl  =  concrete  type  with  entries  either  1,  2, 

or  3,  in  which  1  is  normal  weight  concrete 

(carbonate  aggr.), 

2  is  normal  weight  concrete  (silicon  aggr.), 

and  3  is  lightweight  concrete 
NR2  =  niunber  of  different  categories  of 
^  steel  tendons 

RD3  =  capacity  reduction  factor 

RD4  =  specified  con5>ressive  strength  of 

concrete  (psi) 
RD5  =  emissivity 
RD6  =  convection  heat  transfer  coefficient 

(btu/in2.min.F) .   If  RD5,  RD6  =  0, 

default  values  will  be  used. 

SECTION  FOUR:  BEGINNING  OF  THERMAL  MODEL  DATA 
Line  #1 :  Thermal  control  separator 

THERMAL 
Line  #2 :  Control  data  for  thermal  model 

NRl ,  NR2 ,  NR3 ,  NR4 ,  NR5 ,  NR6 ,  NR7 ,  NR8 

where , 

NRl  =  number  of  nodes  for  the  thermal  model 

NR2  =  number  of  time  steps 

NR3  =  number  of  sides  exposed  to  fire 

NR4  =  number  of  elements  of  the  thermal  model 

NR5  =  mode  of  conductivity  data  input 

NR6  =  mode  of  thermal  capacity  data  input 

NR7  =  number  of  points  defining  conductivity  data 

NR8  =  niimber  of  points  defining  capacity  data 

NRl  =  0  default  values  used, 

=  1  user  input  data,  for  1=5  and  6. 

SECTION  FIVE:  NODAL  COORDINATES  FOR  THERMAL  MODEL 

Line  #1 :  Node  control  header 

NODES 
Subsequent  lines :  nodal  coordinates 

NI  Z=XI  T=TI  G=N1,N2,N3 
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where, 

f  NI  =  node  number 

XI  =  X-  coordinate  for  node  NI 
YI  =  y-  coordinate  for  node  NI 
G  =  initiates  the  generation  sequence 
Nl  =  first  node  in  the  generation  sequence 
N2  =  last  node  in  the  generation  sequence 
N3  =  node  increment  nximber  (default  is  1) 
NOTE:  *  coordinates  of  nodes  Nl  and  N2  must 
have  been  defined 

*  G  is  omitted  if  no  generation 

*  this  section  must  end  with  a  blank  line 

SECTION  SIX:  TIME -TEMPERATURE  RELATIONSHIP 

r^*  ,-      Line  #1 :  Heading  separator 

^  TIMESTEPS 

Subsequent  lines  (NR2  in  Sect.  Four)  :  Time-tertp  data 

NR1,RD2,RD3 

v.!;;  where,    ,   '^    '■  ^• 

NRl  =  time  step  number 

RD2  =  time  (minutes) 

RD3  =  temperature  (degr.  F) 

SECTION  SEVEN:  USER  INPUT  CONDUCTIVITY  DATA  (omit  this 
section  if  NR5  =  0  in  Section  four 

Line  #1 :  Control  header 

CONDUCTIVITY 
Subsequent  lines  (NR7  in  Section  Four) 

RD1,RD2 

where , 

RDl  =  tenperature  value 
RD2  =  conductivity  value 

SECTION  EIGHT:  USER  INPUT  CAPACITY  DATA  (omit  this  section 
if  NR6  =  0  in  Section  four 

Line  #1 :  Control  header 

CAPACITY 
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Subsequent  lines  (NR8  In  Section  Four) 

RD1,RD2 

where , 

RDl  =  tenperature  value 

RD2  =  thermal  capacity  value 

SECTION  NINE:  BOUNDARY  (EXPOSED)  SURFACES  DEFINITION 

Line  #1:  Surface  control  header 

SURFACE 

Subsequent  lines  (NR3  in  Sect.  Four) :  Exposed 
surface  data 

NR1,NR2,NR3 

where , 

NRl  =  side  nuinber 

NR2  =  node  I  of  side  NRl 

NR3  =  node  J  of  side  NRl 

SECTION  TEN:  ELEMENT  INFORMATION 

Line  #1 :  Element  header  separator 

Vl.VMVMT 

Subsequent  lines  (NR4  in  Sect.  Four) :  Element 
description  data 

NRl ,  NR2 ,  NR3 ,  NR4 ,  NR5 ,  NR6 

where , 

NRl  =  element  number 

NR2  =  node  I  of  element  NRl 

NR3  =  node  J 

NR4  =  node  K 

NR5  =  nximber  of  element  sides  exposed  to  fire 

NR6  =  number  of  the  side  parallel  to  x-axis,  where 

1  is  for  side  bound  by  nodes  I-J 

2  is  for  side  bound  by  nodes  J-K 

3  is  for  side  bound  by  nodes  K-I 

NOTE:  IF  ONLY  THERMAL  ANALYSIS  IS  REQUIRED,  (i.e.  NR11=0  IN 
SECT.  ONE),  THE  INPUT  DATA  FILE  MAY  ENDS  HERE. 
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j  ''J  SECTION  ELEVEN:  BEGIN  OF  STRUCTURE  MODEL  DATA 

\     :  .  Line  #1 :  Control  heading  separator 

STRUCTURE 

Line  #2 :  Control  data 

NR1,NR2,NR3 

where , 

NRl  =  number  of  nodes  for  the  structure  model 
NR2  =  number  of  beam  elements 
NR3  =  number  of  steel  tendons 

SECTION  TWELVE:  NODE  INFORMATION 

Line  #1 :  Coordinate  separator 

COORDINATES 

Subsequent  lines:  nodal  coordinates  data 

NI  X=XI  Y=TI  Z=ZI  G=N1,N2,N3 

where  the  terms  are  as  defined  in  Section  Five  above 

NOTE:  *  This  section  must  end  with  a  blank  line 

SECTION  THIRTEEN:   BOUNDARY  CONDITIONS 

Line  #1 :  Boundary  control  header 

BOUNDARY 

Subsequent  lines:  Node  boundary  conditions 

N1,N2,N3   DOF=FDX,FDY,FDZ,FRX,FRY,FRZ 

where , 

Nl  =  first  node  in  the  generation  sequence 
N2  =  last  node  in  the  generation  sequence 
N3  =  nodal  increment  nxuhber  (default  is  1) 
FD's,  FR's  represent  the  boundary  condition  status 
for  the  three  nodal  translations  and 
rotations  respectively.  Each  entry  takes  the 
symbol  either  F  or  R,  in  which, 

R  =  free  translation/  rotation  condition 
F  =  fixed  condition 
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NOTE:  *  If  no  node  generation,  N2  and  N3  are  omitted 
*  This  section  must  end  with  a  blank  line 

SECTION  FOURTEEN:  BEAM  ELEMENT  DATA 

Line  #1 :  Heading  separator 

BEAM 

Subsequent  lines  (NR2  in  Sect.  Eleven) :  Element 
definition 

NR1,NR2,NR3   W=BT    I=HI,BI   J=HJ,BJ 

where , 

NRl  =  beam  element  number 

NR2  =  node  I  of  beam  NRl 

NR3  =  node  J  of  beam  NRl 

BY  =  element  load  (uniformly  distributed  in 

lb/in)  in  the  positive  y-direction 
HX  =  overall  height  of  section  at  end  X  (in.) 
BX  =  width  of  compression  zone  at  end  X  (in.) 

of  beam  element  NRl 

SECTION  FIFTEEN:   INFORMATION  FOR  PRESTRESSING  STEEL 

Line  #1:  Steel  control  header 

STEEL 

Subsequent  lines  (NR3  in  Sect.  Nine) :  Steel 
property  data 

NRl   A=APS   S=FPU   P=PE 

where , 

NRl  =  steel  property  type  number 

APS  =  area  of  steel  tendon  type  NRl 

FPU  =  ultimate  tensile  strength  of  steel 

PE  =  effective  prestressing  force  (after  loses) 

Subsequent  lines  (NR2*NR3  in  Sect.  Nine) :  Segment 
definition 

NR1,NR2,NR3,RD4,RD5   M=MI 

where , 

NRl  =  beam  element  number  containing  the 
steel  segment 


';^i  ;:,  NR2  =  steel  segment  nmnber  in  eleinent  NRl 

'*  *'"^-^  NR3  =  subelement  ntunber  defining  steel  tenperature 

RD4  =  y-eccentricity  of  end  I  of  the  steel  segment 
f^       "  RD5  =  y-eccentricity  of  end  J  of  segment  NR2 

f"  •■  MI  =  material  property  number  of  steel 

segment  NR2 

'  '   SECTION  SIXTEEN:  EXTERNAL  CONCENTRATED  LOADS 

Line  #1:  Loads  control  header 

LOADS 
Subsequent  lines:  Nodal  load  data 

N1,N2,N3   L=LC   F=FX,FY,FZ,FXX,FTY,FZZ 

where , 

Nl ,  N2 ,  N3  bare  the  same  definitions  as 

in  Sect.  Eleven  '-"4 

LC  =  load  case  number 
FX  =  nodal  load  in  x-direction 
FY  =  nodal  load  in  y-direction 
FZ  =  nodal  load  in  z-direction 
FXX  =  moment  about  x-axis 
FYY  =  moment  about  y-axis 
FZZ  =  moment  about  z-axis 

NOTE:  *  If  no  concentrated  loads,  this  section  j 

can  be  omitted  ] 

*  This  section  must  end  with  a  blank  line.  ' 


**  END  ** 
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EXAMPLE  OF  INPUT  DATA  FILE 
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NWBEAM3 
1,1,1,1,1,1,1,1,1,1,1,1 

C: allowable  errors 

CONVERGENCE 
20,20,0.05,0.005,0.01,0.01,20 


■l\  r-  '  1   C: material  input  data 

MATERIALS 
.   '-.   2,2,1.0,5000,0.5,5.09E-04 

\      C:input  data  for  thermal  model 

THERMAL  MODEL  CONTROL  DATA 
116,43,28,176,1,1,4,9 

C : node  coords  for  thermal  model 


NODES 

1  X=0     Y=0 

7  X=6  G=l , 7 

8  X=6 . 75 

9  X=0    Y=0 . 5 

15  X=6  G=9,15 

16  X=6.75 

17  X=0    Y=1.5 

23  X=6  G=17,23 

24  X=6.75 

25  X=0    Y=2.5 

31  X=6  G=25,31 

32  X=6.75 

33  X=0    Y=3 . 5 

39  X=6  G=33,39 

40  X=6.75 

41  X=0    Y=4.25 

47  X=6  G=41,47 

48  X=0    Y=5.25 

53  X=5  G=48,53 

54  X=0     Y=6.25 

58  X=4  G=54,58 
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59  X=0    Y=7.25 

62  X=3 

G=59,62 

63  X=0    y=8 . 25 

78         y=13.25 

G=63,78,3 

64  X=l    Y=8.25 

79         y=13.25 

G=64,79,3 

65  X=2    Y=8 . 25 

80         Y=13.25 

G=65,80,3 

81  X=0    Y=14 

83  X=2 

G=81,83 

84  X=0    Y=15 

87  X=3 

G=84,87 

88  X=0    Y=16 

92  X=4 

G=88,92 

93  X=0    Y=17 

98  X=5 

G=93,98 

99  X=0    Y=18 

104  X=5 

G=99,104 

105  X=0    Y=19 

110  X=5 

G=105,110 

111  X=0    Y=20 

116  X=5 

G=lll,116 

Crtime  temperature 

data 

TIMESTEP   TEMPERATURE  DATA 

1,0,68 

2,5,1000 

3,10,1300 

4,15,1399 

5,20,1462 

6,25,1510 

7,30,1550 

8,35,1584 

9,40,1613 

10,45,1638 

11,50,1661 

12,55,1681 

13,60,1700 

14,65,1718 

15,70,1735 

16,75,1750 

17,80,1765 

18,85,1779 

19,90,1792 

20,95,1804 

21,100,1815 

22,105,1826 

23,110,1835 

24,115,1843 

25,120,1850 

26,130,1862 

•■  ^ 
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27,140,1875 
;  , T  28,150,1888 
29,160,1900 
30,170,1912 
31,180,1925 
32,190,1938 
33,200,1950 
34,210,1962 
35,220,1975 
36,230,1988 
37,240,2000 
38,250,2012 
39,260,2025 
40,270,2038 
41,280,2050 
42,290,2062 
43,300,2075 

C:user  input  thermal  conductivity  data 

CONDUCTIVITY 

68,0.001193 

560,0.001057 

1112,0.000904 

1470,0.000804 

C:user  input  theraml  capacity  data 


CAPACITY 

68,0.015534 

390,0. 

,023349 

752,0, 

,023304 

770,0. 

,024424 

788,0. 

,025546 

815,0. 

,027229 

833,0. 

,028351 

932,0. 

034522 

1112,0.023303 

C: boundary  su 

SURFACE 

1,1,2 

2,2,3 

3,3,4 

4,4,5 

5,5,6 

6,6,7 

7,7,8 

8,8,16 

9,16,24 

10,24, 

,32 
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11,32,40 
12,40,47 
.  13,47,53 
14,53,58 
15,58,62 
16,62,65 

•1  ^  17,65,68 
18,68,71 
19,71,74 

\\'-:  20,74,77 
21,77,80 
22,80,83 
23,83,87 
24,87,92 
25,92,98 
26,98,104 
27,104,110 
28,110,116 


C: element  information 


ELEMENTS 

1,1,2,9,1,1 
^   2,2,10,9,0,2 
fa    3,2,3,10,1,1 
^  ^     4,3,11,10,0,2 

5,3,4,11,1,1 
i^'      6,4,12,11,0,2 

7,4,5,12,1,1 


8,5,13,12,0,2  J 

9,5,6,13,1,1 

10,6,14,13,0,2 

11,6,7,14,1,1 

12,7,15,14,0,2 

13,7,8,15,1,1 

14,8,16,15,1,2 

15,9,10,17,0,1 

16,10,18,17,0,2 

17,10,11,18,0,1 

18,11,19,18,0,2 

19,11,12,19,0,1 

20,12,20,19,0,2 

21,12,13,20,0,1 

22,13,21,20,0,2 

23,13,14,21,0,1 

24,14,22,21,0,2 

25,14,15,22,0,1 

26,15,23,22,0,2 

27,15,16,23,0,1 

28,16,24,23,1,2 

29,17,18,25,0,1  ; 

30,18,26,25,0,2 
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31 

rl8 

,19 

,26 

rO, 

32 

rl9 

,27 

,26 

,0, 

33 

rl9 

,20 

,27 

0, 

34 

,20 

,28 

27 

,0, 

35 

20 

21 

28 

0, 

36 

,21 

29 

28 

0, 

37 

21 

22 

29 

0, 

38 

22 

30 

29 

0, 

39 

22 

23 

30 

0, 

40 

23 

31 

30 

0, 

41 

23 

24 

31 

0, 

42 

24 

32 

31 

1, 

43 

25 

,26 

33 

0, 

44 

26 

34 

33 

0, 

45 

26 

27 

34 

0, 

46 

27 

35 

34 

0, 

47 

27 

28 

35 

0, 

48 

28 

36 

35 

0, 

49 

28 

29 

36 

0, 

50 

29 

37 

36 

0, 

51 

29 

30 

37 

0, 

52 

30 

38 

37 

0, 

53 

30 

31 

38 

0, 

54 

31 

39 

38 

0, 

55 

31 

32 

39 

0, 

56 

32 

40 

39 

1, 

57 

33 

34 

41 

0, 

58 

34 

42 

41, 

0, 

59 

34, 

35 

42 

0, 

60 

35, 

43, 

44, 

0, 

61 

35, 

36, 

43, 

0, 

62 

36 

44, 

43 

0, 

63 

36 

37 

44 

0, 

64 

37 

45, 

44 

0, 

65 

37 

38 

45 

0, 

66 

38, 

46 

45 

0, 

67 

38, 

39 

46 

0, 

68, 

39 

47 

46, 

0, 

69 

40, 

47 

39, 

1, 

70 

41, 

42, 

48 

0, 

71 

42 

49 

48 

0, 

72 

42 

43 

49 

0, 

73 

43 

50 

49 

0, 

74 

43 

44 

50 

0, 

75 

44 

51 

50 

0, 

76 

44 

45 

51 

0, 

77 

45 

52 

51 

0, 

78 

45 

46 

52 

0, 

79 

46 

53 

52 

0, 

80 

47 

53 

46 

1, 

81 

48 

49 

54 

0, 

82 

,49 

55 

54 

0, 

,1 

,2 
,1 
,2 
,1 
,2 
,1 
,2 
,1 
,2 
,1 
,2 
,1 
,2 
,1 
,2 
,1 
,2 
,1 
,2 
,1 
,2 
,1 
,2 
,1 
,2 
,1 
,2 
,1 
,2 
,1 
,2 
,1 
,2 
,1 
,2 
,1 
,2 
/3 
,1 
,2 
,1 
,2 
,1 
,2 
,1 
,2 
,1 
,2 
,3 
,1 
,2 
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!■■#--' 


.^^- 


83,49,50,55, 

0,1 

84,50,56,55,0,2 

85,50,51,56,0,1 

86,51,57,56,0,2 

87,51,52,57,0,1 

88,52,58,57,0,2 

89,53,58,52,1,3 

90,54,55,59,0,1 

91,55,60,59,0,2 

92,55,56,60,0,1 

93,56,61,60,0,2 

94,56,57,61,0,1 

95,57,62,61,0,2 

96,58,62,57,1,3 

97,59,60,63,0,1 

98,60,64,63,0,2 

99,60,61,64,0,1 

100,61,65,64,0,2 

101,62,65,61 

,1,3 

102,63,64,66 

,0,1 

103,64,67,66 

,0,2 

104,64,65,67 

,0,1 

105,65,68,67 

,1,2 

106,66,67,69 

,0,1 

107,67,70,69 

,0,2 

108,67,68,70 

,0,1 

109,68,71,70 

,1,2 

110,69,70,72 

,0,1 

111,70,73,72 

,0,2 

112,70,71,73 

rO,l 

113,71,74,73 

rl,2 

114,72,73,75 

rO,l 

115,73,76,75 

0,2 

116,73,74,76 

0,1 

117,74,77,76 

rl,2 

118,75,76,78 

0,1 

119,76,79,78 

0,2 

120,76,77,79, 

0,1 

121,77,80,79, 

1,2 

122,78,79,81, 

0,1 

123,79,82,81, 

0,2 

124,79,80,82, 

0,1 

125,80,83,82, 

1,2 

126,81,85,84, 

0,2 

127,81,82,85, 

0,1 

128,82,86,85, 

0,2 

129,82,83,86, 

0,1 

130,83,87,86, 

1,2 

131,84,89,88, 

0,2 

132,84,85,89, 

0,1 

133,85,90,89, 

0,2 

134,85,86,90, 

0,1 

^*'- 


135,86,91,90,0,2 

136,86,87,91,0,1 

137,87,92,91,1,2 

138,88,94,93,0,2 

139,88,89,94,0,1 

140,89,95,94,0,2 

141,89,90,95,0,1 

142,90,96,95,0,2 

143,90,91,96,0,1 

144,91,97,96,0,2 

145,91,92,97,0,1 

146,92,98,97,1,2 

147,93,100,99,0,2 

148,93,94,100,0,1 

149,94,101,100,0,2 

150,94,95,101,0,1 

151,95,102,101,0,2 

152,95,96,102,0,1 

153,96,103,102,0,2 

154,96,97,103,0,1 

155,97,104,103,0,2 

156,98,104,97,1,3 

157,99,106,105,0,2 

158,99,100,106,0,1 

159,100,107,106,0,2 

160,100,101,107,0,1 

161,101,108,107,0,2 

162,101,102,108,0,1 

163,102,109,108,0,2 

164,102,103,109,0,1 

165,103,110,109,0,2 

166,104,110,103,1,3 

167,105,112,111,0,2 

168,105,106,112,0,1 

169,106,113,112,0,2 

170,106,107,113,0,1 

171,107,114,113,0,2 

172,107,108,114,0,1 

173,108,115,114,0,2 

174,108,109,115,0,1 

175,109,116,115,0,2 

176,110,116,109,1,3 

C: input  data  for  structure  model 

STRUCTURE 
12,11,5 
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C:node  coords  for  structure  model 

COORDINATES 

1    X=0  Y=0  Z=0 


^  '■:?  ■' :;  ,^- 
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12  X=ll*12        G=l,12 
C:node  boundary  conditions 

BOUNDARY 

I  DOF=F,R,F,F,F,F 
2,10    DOF=R,R,F,F,F,R 

II  DOF=R,F,F,F,F,R 
12      DOF=R,R,F,F,F,R 

C:beain  element  information 


BEAM 

1,1,2 

W=-14 

.926 

1=20.0,10.0 

J=20. 0,10.0 

2,2,3 

W=-14 

.926 

1=20.0,10.0 

J=20. 0,10.0 

3,3,4 

W=-14 

.926 

1=20.0,10.0 

J=20. 0,10.0 

4,4,5 

W=-14 

.926 

1=20.0,10.0 

J=20. 0,10.0 

5,5,6 

W=-14, 

.926 

1=20.0,10.0 

J=20. 0,10.0 

6,6,7 

W=-14. 

.926 

1=20.0,10.0 

J=20. 0,10.0 

7,7,8 

W=-14, 

.926 

1=20.0,10.0 

J=20. 0,10.0 

8,8,9 

W=-14. 

,926 

1=20.0,10.0 

J=20. 0,10.0 

9,9,10 

W=-14, 

,926 

1=20.0,10.0 

J=20. 0,10.0 

10,10,11 

W=-14. 

926 

1=20.0,10.0 

J=20. 0,10.0 

11,11,12 

W=-14. 

926 

1=20.0,10.0 

J=20. 0,10.0 

Crprestressing  steel  information 

STEEL 

1  A=0.080   S=250000   P=11480 

2  A=0.0302  S=250000  P=4305 
1,1,50,3.27,3.27  M=l 
1,2,81,5.74,5.74  M=l 
1,3,50,3.27,3.27  M=l 
1,4,44,3.27,3.27  M=l 
1,5,70,4.55,4.55  M=2 
2,1,50,3.27,3.27  M=l 
2,2,81,5.74,5.74  M=l 
2,3,50,3.27,3.27  M=l 
2,4,44,3.27,3.27  M=l 
2,5,70,4.55,4.55  M=2 
3,1,50,3.27,3.27  M=l 
3,2,81,5.74,5.74  M=l 
3,3,50,3.27,3.27  M=l 
3,4,44,3.27,3.27  M=l 
3,5,70,4.55,4.55  M=2 
4,1,50,3.27,3.27  M=l 
4,2,81,5.74,5.74  M=l 
4,3,50,3.27,3.27  M=l 
4,4,44,3.27,3.27  M=l 
4,5,70,4.55,4.55  M=2 
5,1,50,3.27,3.27  M=l 
5,2,81,5.74,5.74  M=l 
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5 

5 

5 

6 

6 

6 

6 

6 

7 

7 

7 

7 

7 

8 

8 

8 

8 

8 

9 

9 

9 

9 

9 

10,1,50 

10,2,81 

10,3,50 

10,4,44 

10,5,70 

11,1,50 

11,2,81 

11,3,50 

11,4,44 

11,5,70 


,3,50 

,3.27, 

,4,44 

,3.27, 

,5,70 

,4.55, 

,1,50 

,3.27, 

,2,81 

,5.74, 

,3,50 

,3.27, 

,4,44 

,3.27, 

,5,70 

,4.55, 

,1,50 

,3.27, 

,2,81 

,5.74, 

,3,50 

3.27, 

,4,44 

3.27, 

,5,70, 

4.55, 

,1,50, 

3.27, 

,2,81, 

5.74, 

,3,50, 

3.27, 

,4,44, 

3.27, 

,5,70, 

4.55, 

,1,50, 

3.27, 

,2,81, 

5.74, 

,3,50, 

3.27, 

,4,44, 

3.27, 

,5,70, 

4.55, 

27 
27 
55 
27 
74 


3.27 
3.27 


55 
27 
74 
27 
27 
55 


3.27 
5.74 


,3.27 
,5.74 
,3.27 
,3.27 
,4.55 
,3.27 
,5.74 
,3.27 
,3.27 
,4.55 


27 
27 
55 
27 
74 
27 
27 
55 
,3.27 
,5.74 
,3.27 
,3.27 
,4.55 
,3.27 
,5.74 
,3.27 
,3.27 
,4.55 


M=l 
M=l 
M=2 
M=l 
M=l 
M=l 
M=l 
M=2 
M=l 
M=l 
M=l 
M=l 
M=2 
M=l 
M=l 
M=l 
M=l 
M=2 
M=l 
M=l 
M=l 
M=l 
M=2 
M=l 
M=l 
M=l 
M=l 
M=2 
M=l 
M=l 
M=l 
M=l 
M=2 


C: external  point  loads 


LOADS 
1   L=l 
5 
9 


F=0. 0,-1250. 0,0, 0,0,0 
F=0. 0,-2500. 0,0, 0,0,0 
F=0. 0,-2500. 0,0, 0,0,0 


C:end  of  input  data  file 


If.  f 


•r. 
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